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Abstract. We consider compressible pressurelcss fluid flows in Lagrangian coordinates in one space dimen- 
sion. We assume that the fluid self-interacts through a force field generated by the fluid itself. We explain 
how this flow can be described by a differential inclusion on the space of transport maps, in particular when 
a sticky particle dynamics is assumed. We study a discrete particle approximation and we prove global 
existence and stability results for solutions of this system. In the particular case of the Euler-Poisson system 
in the attractive regime our approach yields an explicit representation formula for the solutions. 
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1. Introduction 

In this paper, we consider a simple model for one-dimensional compressible fluid flows under the influence 
of a force field that is generated by the fluid itself. It takes the form of a hyperbolic conservation law for 
the density g, which is a nonncgative measure in time and space and describes the distribution of mass or 
electric charge, and the real- valued Eulcrian velocity field v. For suitable initial data (g, v)(t = 0, •) =: (g, v), 
the unknowns (g,v) satisfy 

d t g + d x (gv) =0 1 
d t {Qv) + d x {gv 2 ) = f[g]\ 



[0,co)xK. (1.1) 



The first equation in (|1.1[) . called the continuity equation, describes the local conservation of mass or 
electric charge. Without loss of generality, we will assume in the following that the total mass/charge is 
equal to one initially and that the quadratic moment is finite so that g(t, •) € ^(R) for all t ^ 0, with ^(R) 
the space of probability measures with finite quadratic moment. The second equation in (|1.1[) describes the 
conservation of momentum. We will assume in the following that v(t, •) G J?? 2 (R, g(t, •)) for alH ^ so the 
kinetic energy is finite. 

The continuous map /: .^(R) — > (K) in (jl.lj) describes the force field, with ^(ffi) the space of all 
signed Borel measures with finite total variation. The force depends on the distribution of mass or electric 
charge and we will assume that f[g] is absolutely continuous with respect to g. For further assumptions see 
Section [SJ The typical (simplest) form of / is 

f[o\ = -gd x q g with q e (x) = V(x) + fw(x- y) dg(y) (1.2) 

for suitable C 1 potential functions V, W with (at most) linearly growing derivatives. 
Another relevant example we have in mind is the Euler-Poisson system, for which 

f[o\ = ~Qd x q e with q g solution of — d xx q e = X(g — a). (1-3) 

When g is absolutely continuous with respect to the one-dimensional Lebesgue measure C 1 , then the function 
q e admits a representation similar to (|1.2p . with 

V(x):=~ f\x-y\da(y), W(x) := ±\x\. (1.4) 
L Jr. z 

If p is not absolutely continuous with respect to C , then we have a similar representation with a nondiffcr- 
entiable W, so that f[p] must be defined by a suitable approximation process. 

The Euler-Poisson equations in the repulsive regime (with A < and negative concave potential W) is a 
simple model for semiconductors. In this case, g describes the electron or hole distribution and the scalar 
function q e represents the electric potential generated by the distribution of charges in the material. Here a 
is the concentration of ionized impurities. The Euler-Poisson equations in the attractive regime (with A > 
and positive convex potential W) is the one-dimensional version of a cosmological model for the universe 
at an early stage, describing the formation of galaxies. Now q g represents the gravitational potential and 

(7 = 0. 



1.1. Singular solutions and particle models. Since there is no pressure in (jl.ll) . there is no mechanism 
that forces the density g to be absolutely continuous with respect to the Lebesgue measure. In fact, the 
system p. II) admits solutions that are singular measures. Assume that we are given initial data in the form 
of a finite linear combination of Dirac measures: 

N N 

g = *^2fhiS Xi and gv = ^ fhiViS Xi , (1.5) 

i=l i=l 

where x = (xi, • • • , xn) £ K w are the initial locations of N particles denoted Pi, ■ ■ ■ , Pn, with corresponding 
masses fh = (mi, • • • , m_/v) and initial velocities v = (v\, ■ ■ ■ , Vn). We require that fhi > and J2i = 1 
so that g € ^(R). For all times t ^ 0, we can assume that the positions x(t) = (xi (£),■■• , xjv(i)) arc 
monotonically ordered, so that they are unambiguously determined and attached to the particles. Then (at 
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least formally) there is a solution of p. II) in the form of a linear combination of Dirac measures: 

N N 

•) = X) m »*W) and (HCM = ^2miVi(t)S Xi ( t ), (1.6) 

i=l i=l 
where the functions (jCi, Uj) solve the system of ordinary differential equations 

±i(t) = Vi(t), Vi(t) = Ojn,i(x(t)) and (xi,Vi)(t = 0) = (xi,Vi) (1.7) 

between particle collisions. Here a m ,i(&) is the value in the point Xi(t) of the Radon-Nikodym derivative of 
the force f[g] with respect to the measure g, so that 

N N 

f[g] =^a mti (x)mi6 Xi if g = ^ m i &xi > C 1 - 8 ) 

i=l i=l 

which is well defined when all N particles are distinct. 

Upon collision of, say, two particles with masses m k and mk+i at some time t > 0, the velocities of each 
one of them are changed to 

f , , x ,, , x m k v k (t-) + m k+1 v k+1 (t-) 

v k {t+) = v k+1 (t+) = ■ , (1.9) 

m k + m k+ i 

so that the momentum is preserved during the collision. Since both particles continue their journey with the 
same velocity, they may be considered as one bigger particle with mass m k + m k +i- Collisions of more than 
two particles can be handled in a similar fashion. We will refer to any solution of (|1.1[) in the form (|1.6j) as 
a discrete particle solution and we will say that it satisfies a global sticky condition if particles after collision 
are not allowed to split. In this case, after each collision, one could relabel the particles so that the system 
(|1.7[) still makes sense (with N reduced in each particle collision) and induces a global in time evolution. 
Let us denote by K N the closed cone 

K N := {x £ R N : x x s£ x 2 x N }, (1.10) 

whose interior is intK^ = {a; 6 M : X\ < x-i < ■ ■ ■ < x^}- The construction of discrete particle solutions as 
outlined above can be done rigorously whenever the functions a m ^ : hit K w — > M. N are uniformly continuous 
in each bounded set (so that they admit a continuous extension to K w still denoted by a m: i) and satisfies 
the compatibility condition 

a m ,k( x ) = a m-k+ i(x) if x k = x k+1 for some 1 < k < N. (1-11) 

This is certainly the case when the potentials V, W considered in (|1.2[) are of class C 1 . On the other hand, 
the case of the Eulcr-Poisson system is much more subtle and presents different features in the attractive or 
the repulsive case. 

The Euler-Poisson case in the repulsive regime: splitting and collapsing of masses. Let us consider the 
simplest situation of N distinct particles with equal initial velocities v, in the repulsive regime A = — 1 with 
(7 = 0. Let M := 0, M t := Y!j=i m j for i = 1, . . . , JV and set A; := §(Mj_i + M, - 1). Then it is not 
difficult to check (see Example I6.9[) that in the repulsive regime 

a mi i(x) = Ai for alH if x G xntK N , (1.12) 

and so there is no continuous extension satisfying p. lip . Starting from distinct initial positions, particles 
follow (at least for a small time interval) the free motion paths 

Xi(t)=Xi + tv + ^Ait 2 . (1.13) 

Since Ai ^ Aj+i for all i, there are no collisions. Taking the limit as the initial positions of two or more 
particles coincide we obtain the same representation for every x € K. N . On the other hand, if two particles 
P k ,P k+ i coincide at the time t = 0, i.e. x k = x k +\ = x with the same initial velocity v, then the "sticky" 
solution x k (t) = x k+ i(t) = x + tv + \{A k + A k+ i)t 2 gives raise to an admissible solution to (jl.ip which is 
different from the previous one. One could also consider a solution where P k ,P k +\ stick in a small initial 
time interval [0, s] and then evolve according to (|1.13|) . 

Considering a situation where the number TV of admissible particles grows to infinity with a uniform initial 
mass distribution concentrating at the origin, one can guess that a "repulsive" solution arising from a unit 
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mass concentrated at x should istantaneously diffuse, becoming absolutely continuous with respect to the 
Lebesgue measure C 1 : the explicit formula is 

g(t,-) =u(t,-)£ 1 with u(t,x) = -^X( S +vt_fc2 iS+ijt+ fr2)(x) for all x 6 R and t ^ 0. (1.14) 

An even more complicated situation occurs e.g. if w,; = for i =/= k, k + 1, but Vk > > Vk+i in such a way 
that a collision occurs between Pk and Pk+i at some time t = r, after which the particles could stick or wait 
for some time and then evolve as in the previous example. 

It would be important to find a selection mechanism that gives raise to a stable notion of solution and 
to obtain a continuous model by passing to the limit in the number of particles. In this paper we study a 
criterium of the following type : Assume that two particle Pk,Pk+i collide at some time r > with incoming 
velocities t>fc(r_) ^ Vk+i(r~). Then the particles will stick together for all times r < t < s provided that s 
is small enough so that 

Vk(r-)+ a m . k (x(T))dT > v k (t) = v k +i(t) ^ v k +i{r-) + / a m ,fe + i(x(r)) dr (1-15) 

J r J r 

for all r < t < s. Conversely, if (|1.15|) becomes false for some time s > r, then the particles may separate 
again. A rigorous formulation of condition (|1.15[) in the case of a simultaneous collision or separation of more 
than two particles, or of a continuous distribution of masses, can be better understood in the famework of 
differential inclusions in a Lagrangian setting, which we will describe in Section |5. II Before giving an idea 
of this approach, let us brefly consider how (| 1 . 1 5[) greatly simplifies in the attractive regime. 

The attractive Euler-Poisson system and the sticky condition. In the attractive case, we can simply invert 
the signs in (|1.12[) . It turns out, however, that the behaviour of the two-particles example considered in the 
previous paragraph changes completely, since the limit when two particles collapse exhibit a strong stability: 
after a collision, two or more particles stick together and do not split ever again, giving raise to a global 
sticky solution. 

This reflects the fact that the sticky condition in the attractive regime implies (| 1 . 1 5[) for all s > r: the 
functions a m ^ defined by the negative of (|1.12p always satisfy a m ^k{x) ^ a m _k+i{x) and the incoming 
velocities of two particles Pk,Pk+i colliding at some time r satisfies i>fc(r— ) ffe+i (?*-), so that any sticky 
evolution corresponding to Xk{t) = Xk+i(t) for t ^ r will satisfy (|1.15[) . 

As we shall see, the differential description in the Lagrangian setting we will adopt encodes (|1.15|) and 
corresponds to a sticky condition whenever the acceleration field is continuous (as in (|1.11[0 or it is of 
attractive type. In the repulsive case it will model a suitable relaxation mechanism allowing for separation 
of particles after collision, still preserving the stability of the evolution. 

1.2. Lagrangian description and differential inclusions. In this paper, we will give an interpretation 
of system (|1.1[) m the framework of differential inclusions. As before, let us first consider the simpler case 
of the dynamic of a finite number of particles. We can identify the positions of a collection of particles 
Pi, ■ ■ ■ , Pjv with a vector x = (xi, • • ■ , xn) € M. N : since we labeled the particles in a monotone way, it is 
not admissible for particles to pass by one another, so the order of the locations must be preserved and the 
vector x is confined in the closed convex cone K N defined by (|1.10[) . Denoting by v = (vi, ■ ■ ■ , Vn) £ M. N 
the vector of the velocities of the particles, their trajectories between collisions are mostly determined by a 
system of differential equation 

x(t)=v{t), v{t) = a m (x(t)), (1.16) 

where a m {x) := (a r7li i(x), • • • ,a m ^(x)) is a vector field defined for x S K N as in (|1.7p . which in the 
simplest case is continuous. Whenever the vector x(t) hits the boundary of the domain 

dK N :={xeK N -M x ^0}, Sl m :={j:xj=x j+1 ,j = l...N-l}, (1.17) 

however, an instantaneous force changes its velocity in such a way that it stays inside of K N . 

In order to find a mathematical model that describes this situation, we must first identify the set of 
admissible velocities at each point x £ K N , which is called the tangent cone of K N at x. It is defined by 

T X K N := cl{% - x) : y e K N , 9^ oj. (1.18) 
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Figure 1 . Projection of velocities onto the tangent cone. 



In our situation, it is not difficult to check that 

T X K N = ER N : Vj sC v j+ i for all j G fl x j. (1.19) 

Identity (|1.19[) shows that when two particles collide, then the velocity of the left particle cannot be greater 
than the velocity of the right particle, so the left particle cannot pass. 

Assume now that x{t) G dK N at some time t and let v(t— ) be the velocity immediately before the impact. 
That is, let v(t—) be the left-derivative of the curve 1 1-> x(t) at time t. The instantaneous force that is active 
on impact must change the velocity to a new value in the tangent cone T X ^K. N of admissible velocities. 
Typically, there are many possibilities. Assuming inelastic collisions, we impose the impact law: 

v(t+) := P Tm(t)K Nv(t-), (1.20) 

where v(t+) is the velocity immediately after impact (the right-derivative of x(t)). We denote by Pt^^k" 
the metric projection onto T x u-\M. so that 

\\v(t-)-v(t+)\\ = min{\\v(t-)-u\\ m :ueT x{t) K N y 

Hence v(t+) is the element in T X ^~K N closest to v(t— ) with respect to the weighted Euclidean distance 
induced by the norm 



N 



mivf for all v G T x(t) K N c R N , (1.21) 

\ i=i 

and therefore unique: see Figure [1] for a graphic representation of this rule. 

It is well-known that the metric projection onto closed convex cones admits a variational characterization 
of its minimizers; see |24| . In particular, we have 



(v(t-) - «(*+)) • u ^ for all u G T x(t) J 



-N 



We deduce that the instantaneous force that changes the velocity upon impact onto the boundary OK 
must be an element of the normal cone N X ^1K N , which is defined as 



/v 



N^K N 



jn G R N : n ■ (y - x) ^ for all y G K w |. (1.22) 

Note that the normal cone N X K N equals the subdiffcrential dl^N (x) of the indicator function I k n of K N at 
the point x. This follows immediately from the definition of the subdifferential. 
This suggests to consider the second-order differential inclusion 

x = v, v + N X K N 3 a m {x) in[0,oo). (1.23) 

Notice that since v can exhibit jumps, solutions to (|1.23[) should be properly defined in a weak sense in the 
framework of functions of bounded variation. Second-order differential inclusion have been studied in the 
literature and existence of solutions has been shown in a genuinely finite dimensional setting. We refer the 
reader to 15, 15, 21^ and the references therein for further information. Due to the possible nonuniqueness of 
solutions to second-order differential inclusions [5T] and to the lack of estimates to pass to the limit when 
N — 5> oo, we need a better understanding of the particular features of our setting, in particular of the convex 
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The sticky condition and an equivalent formulation of (jl.23l) . It has been shown in |17l that the one- 
parameter family of normal cones N X ^K N along an evolution curve x : [0, oo) — > WL N for which a gobal 
stickyness condition holds, satisfies the remarkable monotonicity property 

N x{s) K N C N x{t) K N for all s < t. (1.24) 

Consequently, for any selection £ : [0, oo) — > M. N satisfying £(£) G N X ^K N (such as a m (x(t)) — v{t) in 
([L23]) ) we have 

J £(r) dr G N^K* for all s < f. 
An integration of (| 1 . 23[) yields, at least formally, that 

v(t) + N x(t) K N 3 v(s) + j a(x(r)) dr 

and therefore the system (|1.23j) can be rewritten in the form 

x = v, v + N X K N 3 y, y = a m (x). (1.25) 

Introducing new unknowns (x,y), we can rewrite (|1.25[) as a first order evolution inclusion 

x + N x(t) K N 3 y 

() y (1.26) 
y = a m (x) 

for which an existence and stability theory is available, at least when Lipschitz map. 

We will show that formulation (|1.26|) enjoys interesting features and always induces a measure-valued 
solution to (|1.1[) . When the field a m satisfies the compatibility condition (jl.lll) . solutions to (|1.26l) satisfies 
the sticky condition, and the same property holds also for the Eulcr-Poisson equation in the attractive 
regime. In the repulsive case, we will see that (|1.26[) is a robust formulation of condition fj 1 . 1 5|) . Let us now 
consider the infinite-dimensional case. 



1.3. Diffuse measures and differential inclusions for Lagrangian parametrizations. In order to 
deal with general measure- valued solutions of (|1.1[) , we had to recourse to Lagrangian coordinates, using 
ideas of optimal transport as considered in [17] . 

Monotone Lagrangian rearrangemens. In this approach, the discrete set of parameters {1, 2, • • • , N} involved 
in the representation of discrete particle measures (|1.5[) will be substituted by f2 = (0, 1). For every particle 
labeled by m £ H, we will denote by X(t,m) G K its position at time t. The map X can be uniquely 
characterized in terms of the measure g: it is the uniquely determined nondecreasing and right-continuous 
map X : O — > M such that 

X(m) ^ x m< g((-oo,x}) for all (1.27) 

Equivalently, the push-forward X#m of the one-dimensional Lebesgue measure m := C 1 ^ under the map X 
equals g. Recall that the push-forward measure is defined by 

X # m(A) := miX^iA)) for all Borel sets A c R. (1.28) 

Therefore the map X is the optimal transport map pushing m forward to g. We refer the reader to Section [2] 
for further explanation. 

In this way, to any solution (g,v) of (jl.lj) . we can associate a map X : [0, oo) x f2 — > R with X(t,-) 
nondecreasing and a velocity V : [0, +oo) x SI — s> R such that 

X{t, -)#m = g(t, •) foralU^O, V(t, ■) = v{t, X(t, ■)) = d t X(t, ■). (1.29) 

Our goal is to show that (jl.ip can be associated to a differential inclusion in terms of (X, V). This observation 
allows us to derive existence and stability results (see Sections [3] and H]) for (suitably defined) solutions of 
(|l.ip . which together with the existence of discrete particle solutions (see Section [5]) imply a global existence 
result for (11.11) for general initial data; see Section O 
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Differential inclusions. The framework of first-order differential inclusions, analogous to the setting we al- 
ready discussed for the discrete case (|1.23|) , serves as a guiding principle for our discussion. The role of the 
cone is now played by the cone of optimal transport maps 

X := (X E ^ 2 {VL) : A is nondecreasing j (1.30) 

in the Hilbert space J$? := Jz? 2 (51). Even if in this infinite dimensional setting the boundary of is dense, 
we can still consider the normal cone Nx^ for given X G J(f , which is given by 

N x J(f := |w G Jz? 2 (ft) : J w(x - a) dm < for all X <= jr}. (1.31) 

Again we have that = dlj^(X). It can be shown that = {0} if and only if the map X is not 

strictly increasing in fi. That is, whenever Q\ <ZS where 

fix '■= | w € O: X is constant in a neighborhood of mX. (1-32) 

Note that fix is the complement of the support of the distributional derivative of X . 

Consider now a family of densities t H- g(t, •) that satisfies (jl.ljl . Let t M> A(t, •) G be the associated 
family of optimal transport maps; see (|1.29[) . We want to interpret X as a solution of differential inclusions, 
similar to (fQ3|) and (TQ6| . 

Even at the continuous level, the monotonicity property (|1.24[) for sticky particle evolutions plays a crucial 
role. Note that the optimal transport map X G J(T takes a constant value ielon some interval (a, /3) C O 
if the mass /3 — a (the Lcbcsgue measure of the interval) is moving to the same location, thereby forming a 
Dirac measure at x. Therefore sticky evolutions will be characterized as curves t i— > X(t, •) with the property 
that 

for any t\ ^ t% we have Q, X (ti) C Cl X (t 2 )- (1.33) 

Notice that (|1.33|) implies that once a Dirac measure is formed, it may accrete more mass over time, but it 
can never lose mass. It also implies the following statement: It is not possible for mass to jump from one 
side of a Dirac measure to the opposite side. Whenever mass is crossing a Dirac measure, it gets absorbed. 

A formulation via differential inclusions needs a Lagrangian expression of the force term in fll.ip . That 
is, we must find a map F: Jlf — > «5f 2 (J7) with the property that 

/ if)(x) f[g](dx) = / ip(X{m))F[X](m)dm for all ip € 0(R), (1.34) 

JR JQ 

whenever X G J{f and X#m = g. We refer the reader to Section [5] for further discussion about the existence 
and properties of maps F satisfying f| 1 . 34|) . In the following, we will assume that F is continuous as a map 
of JtT into £? 2 (n). 

We then could expect A to be a solution of a second-order differential inclusion, but arguing as for the 
discrete case (|1.23|) at least in the case of sticky evolutions (|1.33l) we end up with 

X(t) +dlje(X(t)) 3 V+ [ F[X(s)]ds (1.35) 

J o 

for a.c. t > 0. This formulation and its consequences is at the heart of our argument. 

It is a remarkable fact (see Theorem I3.5[) that solutions to ()1.35[) always parametrize measure- valued 
solutions to the partial differential equation (|1.1[) . Provided F satisfies suitable continuity properties, it will 
be possible to prove existence (and uniqueness, when F is Lipschitz) of solutions to (|1.35[) for any initial 
data (A, V) G Jff x Jz? 2 (51) by combining the theory of gradient flows of convex functionals in Hilbert spaces 
[5] with suitable compactness arguments. 

When F satisfies a suitable sticking condition, which is satisfied e.g. in the case of C 1 potentials in (|1.2I) 
and of the Eulcr-Poisson system in the attractive regime, then solutions to (|1.35[) form a semigroup and have 
the sticky evolution property (|1.33[) . Even for general F (and in particular for the Euler-Poisson system 
in the repulsive regime) the differential inclusion (| 1 . 35[) still selects a stable parametrization of solutions to 
(|1.1[) . This is somewhat surprising since the reduction from second-order to first-order differential inclusion 
was motivated by the monotonicity (jl.33|) . which typically is false without additional assumptions on F. In 
this Introduction we refer to such solutions as "robust." 
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Representation formulae for the Euler-Poisson system. In the case of the Eulcr-Poisson system (jl.3[) with 
cr = one can show that the Lagrangian representation of the force / is given by 

F[X](m) = -XA(m) where A(m) := m - -. (1.36) 
Note that the map F[X] is independent of X and (|1.35[) becomes 

X(t) + dIj(r{X{t)) BV-XtA 

In the attractive regime when A > 0, an explicit representation formula for the Lagrangian solution can be 
obtained (see Theorems l6.1ip . In fact a careful analysis shows that the solution X to (|1.35[) can be computed 
by solving the trivial ODE in X obtained by eliminating the ^-constraint: 

^-X(t) = V-XtA 
dt v ' 

(whose solution is X{t) = X + t V — ^Xt 2 A), and then projecting X on Jff: 

X(t) = Pjf{X(t)) = Px (x + tV- ^Xt 2 A 

Applying the characterization given in |17j , the metric projection of Pj%r onto can be found by introducing 
the primitive functions 

rrn nm -t 

X(m):= X(£)d£, V(m) := / V(£)d£, A(m) := -{m 2 - m), 

.In .In 2 



and the time evolution 



X(t, m) := X(m) + t V(m) - ^Xt 2 A(m), 



and then taking the derivative with respect to m of the convex envelope X**(t, ■) of X(t, •): 

X(t,m) = ^-X**{t,m), 

which defines a density g t = X(t)#m. It is then a simple exercise to recover formula (|1.14[) in the case 
X(m) = x, V(m) = v, since X(t, m) = X** (t, m) and X(t, m) = x + tv + \t 2 {m — 3). 

1.4. Time discrete schemes. In this section, we show that the first-order differential inclusion (|1.35j) can 
be used to design a stable explicit numerical scheme to compute robust solutions to (|1.1[) . In fact, this 
scheme is essentially the same as the one introduced in [7] for "order-preserving vibrating strings" and 
"sticky particles" , with just mild modifications. For simplicity, we concentrate on the pressureless repulsive 
Euler-Poisson system with a neutralizing background 

d t g + d x (gv) = 0, (1.37) 

dt{gv) + d x {gv 2 ) = -gd x q e - d xx q e = -(g - 1), (1.38) 

(cf. (|1.3[) with A = — 1 and a = 1). We assume the initial conditions to be 1-periodic in x and the density 
g to have unit mean so that the system is globally neutral and the electric potential q e is 1-periodic in x. 
Note that we choose the periodic setting only for convenience. In fact, for any non-periodic solution g of 
(jl.ip , one could consider the push- forward of g under the map x H > x — [x] for all with [x] the largest 

integer not greater than x. One obtains a new density g* that is concentrated on [0, 1) and therefore can be 
extended 1-periodically to the whole real line. One can then show that g* satisfies the same equation. We 
refer the reader to [12] for details. 

For smooth solutions without mass concentration, written in mass coordinates 

V(t, m) = X(t, m) = v(t, X(t, m)), d m X(t, m)g(t, X(t, m)) = 1 

(which requires that d m X(t,m) ^ 0), one can show that the whole system reduces to a collection of inde- 
pendent linear pendulums labeled by their equilibrium position m and subject to 

X(t,m) = V(t,m), V(t, m) + X(t, m) - m = 0. (1.39) 

(Notice that, due to the spatial periodicity of the initial conditions, the new unknown X(t,m) — m and 
V(t, m) arc 1-pcriodic in m.) This reduction is valid as long as the pendulums stay "well-ordered" and do 
not cross each other, i.e., as long as X(t,m) stays monotonically nondecreasing in m. This "non-crossing" 
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condition is not sustainable for large initial conditions and collision generally occur in finite time. To handle 
sticky collisions, the concept of robust solutions introduced in Section [TT31 is a good way to obtain a well-posed 
mathematical model beyond collisions. 

We are now ready to describe the semi-discrete scheme. Given a time step r > and suitable initial 
data (X, V) =: (X t0 , V t ,q), we denote by (X T n (m), V Ttn (m)) the approximate solution at time t n := nr, for 
n = 0, 1, 2, . . ., defined in two steps as follows: 

(1) Predictor step: we first integrate the ODE (|1.39[) and get [/ T ,n+i and X TyTl+ i accordingly 

X T .n+i(m) = to + (X T<n (m) — m) cos(r) + V T , n (m) sin(r), (1-40) 
V Tin +i(m) = -(X Tin (m) - to) sin(r) + V^„(to) cos(r). (1-41) 

(2) Corrector step: we rearrange X Tj „ + i(m) in nondecrcasing order with respect to to and obtain 
X Ttn+ i(m) . Because of the periodic boundary conditions, we have to perform this step with care. 
We rely on the existence, for each map to t-> Y(m) such that Y(m) — to is 1-pcriodic and locally 
Lebcsgue integrablc, of a unique map to n> Y*(m) such that Y*(m) is nondecrcasing in m and 

/ r)(Y*(m)) dm = / r)(Y(m))dm 
Jo Jo 

for all continuous 1-periodic function r\. 

This time discrete scheme becomes a fully discrete scheme, if the initial data X T fi(m) — m and V T ,o(m) are 
piecewise constant on a uniform cartesian grid with step h. (We just have to be careful with the corrector 
step, by using a suitable sorting algorithm for periodic data.) 

To illustrate the scheme, we show the numerical solutions corresponding to initial conditions 

X (m) = to, Vo{m) = 4sin(27TTO). (1.42) 

We use 400 equally spaced grid points m (which corresonds to 400 "well-ordered" pendulums with to as 
equilibrium position) and 5000 time steps (see Figures : 

t = 0.001, 0.005, 0.01, 

so that the final time T of observation is respectively given by 

T = 5, 25, 50. 

On each picture, we show the space-time trajectories of 50 of the 400 pendulums, with space coordinate on the 
horizontal axis and time coordinate on the vertical one. On these pictures, we observe a strong concentration, 
with sticky collisions, of the pendulums at a very early stage (up to time t = 7r/2) around x = 0.5. Later 
on, some pendulums start to unstick and detach from each other (which allows new concentrations at later 
times t 7r around x = and x = 1). Much later, after t = IOtt, there is no further dissipation of energy, 
and, as pendulums touch each other, they always do so with zero relative speed. Then the corrector step is 
no longer active, and the scheme becomes exact (due to the exact integration of the predictor step). At this 
late stage, the solution becomes 27r-pcriodic in time. We study the convergence of the scheme in Section [7] 

1.5. Plan of the paper. We collect in Section[2]a few basic results on optimal transport in one dimension, 
on convex analysis (concerning in particular the properties of the convex cone J€) , and on convex functionals 

In Section [3l after a brief discussion of the basic properties of the Lagrangian force functional F, we 
introduce the notion of Lagrangian solutions to the differential inclusion (|1.35p . Theorem 13.51 collects their 
main properties, in particular in connection with measure- valued solutions to (jl.ip . Sections 13.31 and 13.51 
provide the main existence, uniqueness, and stability results for Lagrangian solutions, whereas Section 13.41 
is devoted to the particular case of sticky evolutions. 

We study in Section 2] a different class of solutions to (|1.35|) , still linked to (jl.ip , that naturally arise as 
limit of sticky particle systems when F does not obey the sticking condition. These solutions exhibit better 
semigroup properties than the Lagrangian solutions introduced in Section [3J but lack uniqueness. 

Section [5] we carefully study the dynamics of discrete particle systems, which we already briefly discussed 
in the Introduction. Discrete Lagrangian solutions associated to systems like (|1.26p are treated in Section lSTl 
where we also show that they can be used to approximate any continuous Lagrangian solution, as the one 
considered in Section [3] The sticky dynamic at the particle level is considered in tj5.2l the main Theorem 
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Figure 2. Space-time trajectories of pendulums, with timestep r = 0.001. 

15.21 provides the basic results, which allow us to replace second-order with first-order evolution inclusion at 
the discrete level and to get sticky evolutions for sticking forces. The particle approach is a crucial step of 
our analysis, since it avoids many technical difficulties arising at the continuous level. The general idea is to 
prove fine properties of the solutions (such as the monotonicity (| 1 . 33[) in the sticking case or a representation 
formula) at the discrete level and then to extend them to the general case by applying suitable stability 
results with respect to the initial conditions. Those are typically obtained by applying contraction estimates 
(in the case when F is Lipschitz) or compactness via Holly's Theorem, by exploiting higher intcgrability and 
monotonicity of transport maps. 

Section [6] applies the Lagrangian formulation to (jl.l[) . presenting some existence and stability results for 
solutions in the Eulerian formalism. 

In Section [7] we prove the convergence of the time discrete scheme introduced in Section 11.41 

2. Preliminaries 

Let us first gather some definitions and results that will be needed later. 

2.1. Optimal Transport. We denote by ^(R m ) the space of all Borel probability measures on R m . The 
push- forward v := Y#[i of a given measure [i G &(M. m ) under a Borel map Y : R m — > R" is the measure 
defined by v(A) := ^l{Y^ 1 {A)) for all Borel sets A C R". We will repeatedly use the change-of- variable 
formula 

/ CO/) (Y#»)(dy) = [ £(Y(x))»(dx), (2.1) 
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Figure 3. Space-time trajectories of pendulums, with timestep r = 0.005. 



which holds for all Borel maps £: R™ — > [0, oo]. 

We denote by .^(R™) the space of all Borel probability measures g £ £P(M. n ) with finite quadratic 
moment: L ra \x\ 2 g(dx) < oo. The Jz? 2 Kantorovich-Rubinstein-Wasserstein distance W2{qi, Q2) between two 
measures Q\,Q2 S .^(R") can be defined in terms of couplings, i.e. of probability measures g £ 3 g (M n x R") 
satisfying 7r#f? = Qi for i = 1 ... 2, by the formula 

W|(ei,ea):=min( / |z - y| 2 (?(dx, dy) : (? e £»(R" x R"), tt^ = A (2.2) 

Here 7^(2:1, X2) := Xi is the projection on the ith coordinate. It can be shown that there always exists an 
optimal transport plan g for which the inf in (|2.2I) is in fact attained. We denote by r opt (ei, Q2) the set of 
optimal transport plans. 

In the one-dimensional case n = 1, there exists a unique coupling g £ r opt (pi, Q2) realizing the minimum 
of (|2.2p (at least when the cost is finite). It can be explicitly characterized by inverting the distribution 
functions of g\ and gi'- for any g £ ^(R) we consider its cumulative distribution function, which is defined 
as 

M e (x) := £>((-oo, x]) for all (2.3) 
Note that then g = d x M e in 5?'(R). Its monotone rearrangement is given by 



X g (m) := inf <^ x: M e (x) > m > for all m £ SI, 



(2.4) 
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Figure 4. Space-time trajectories of pendulums, with timestep r = 0.01. 



where := (0, 1). The map X e is right-continuous and nondecreasing. We have 

m-^C 1 ^, (X e ) # m = g and / ((x) g(dx) = [ ((X,(m)) dm (2.5) 

for all Borel maps (: R — > [0, oo]. In particular, we have that g & ^(R) if and only if X g G «Sf 2 (Sl). The 
Hoeffding-Frechet theorem |T8l Section 3.1] shows that the joint map X gi}B2 : SI — > Ixl defined by 

X gite2 (m) (x ei (to), Xg 2 (m)j for all m E tt, 

characterizes the optimal coupling g £ r op t(gi, Q2) by the formula 

Q= (X eue2 ) # m; (2.6) 
see [H1IH1H3] f° r further information. As a consequence, we obtain that 

Wi( Ql ,Q 2 )= [ \X ei (m)-X g2 (m)\ 2 dm=\\X ei -X g2 \\% Hn) . (2.7) 

The map g *-> X e is an isometry between .^(R) and J?T, where J(f C Jz? 2 (f2) is the set of nondecreasing 
functions. Without loss of generality, we may consider precise representatives of nondecreasing functions 
only, which are defined everywhere. 
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2.2. Some Tools of Convex Analysis for J€ . Let Jfc be the collection of right-continuous nondccrcasing 
functions in Jz? 2 (f2) introduced in (|1.30[) . Then one can check that is a closed convex cone in the Hilbcrt 
space 5? 2 (n). 

Metric Projection and Indicator Function. It is well-known that the metric projection onto a nonempty 
closed convex set of an Hilbert space is a well defined Lipstchitz map (see e.g. [24]): we denote it by 
Pjr : Jz? 2 (ft) — > JT. For all X G £? 2 (fl) it is characterized by 

Y=P^(X) <=► FgJT, \\Y-X\\<? 2m ^\\Y-X\\j? 2(n) for every Y G JT, 

or, equivalently, by the following families of variational inequalities 

Y = P X {X) ^==> Y&Jt, f(X-Y)(Y-Y) < for all Ye X. (2.8) 

Jn 

Pjf(X) admits a more explicit characterization in terms of the convex envelope of the primitive of X 
[T71 Theorem 3.1]: 

Y = P Jt (X) = ^X**(m), X(m)'.= J X(e)d£, (2.9) 
where ^ denotes the right derivative and 

X**(m) := sup |a + bm: a, b 6t, a + 6m<#(m) Vme(0,l)} (2.10) 
is the greatest convex and l.s.c. function below X. 

Let now Ij^ : J>f 2 (5!) — > [0, +oo] be the indicator function of Jf, defined as 



if X G JT, 
+cxd otherwise, 



which is convex and lower semicontinuous. Its subdifferential is given by 

dIjr(X) := jzGJ^ 2 (ft): Ijr(X) >I X {X) + Z(X-X) for all X G „Sf 2 (i7)| (2.11) 

and it is a maximal monotone operator in «if 2 (17); in particular its graph is strongly- weakly closed in 
_Sf 2 (f2) x Jz? 2 (0). Notice that dl, x {X) = for all X J? since in this case I, X {X) = oo; whenever X & Jfr 
we find that 



dIjt(X) = \z G Jz" 2 (0) : ^ ^ Z(X - X) for all X G jr}, 



so that dIj(r(X) coincides with the normal cone defined by (ll.31|) . 

(|2.8p implies the following equivalence: For all X, Y G JS^ 2 (O) we have 



(2.12) 



Y = Pjt(X) ^ X - Y G dIjt(Y). (2.13) 
Decomposing X = Y + Z in (j2~T3|) with Y", Z G Jz? 2 (^), we find that 

zedijtriy) <=► r = p^(r + z). (2.14) 

Lemma 2.1 (Contraction). Lef K — ► (— oo,+cx)] 6e a convex, lower semicontinuous function. For all 
Xi,X 2 G 2 we then have 

r ^(Pjr(^i) - PjrP^)) < / ^(Xi - X 2 ). 

Jo 

In particular, the metric projection Px is a contraction with respect to the I£ v (SX)-norm with p G [1, oo] and 
for allX\,Xi G J2? 2 (f2) we can estimate 

\\Px{X x ) - Pjtr{X 2 )\\x Hn) < - X 2 ||^ P(n) . (2.15) 
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We refer the reader to Theorem 3.1 in [T7] for a proof. Notice that by choosing A 2 = in Lemma [2.11 
for which P,y^(X 2 ) = 0, we obtain the inequalities 

tp(Pjt(X)) ^ [ xl>(X) for all X G ,S? 2 (fi), (2.16) 
Jn 

l|PjrPOII*p(n) ^ ||A||^ (0) for all X G .i? 2 (fi). (2.17) 
A similar result holds for the Jz? 2 -orthogonal projection Pjff x onto the closed subspace Jtx, A G Jt, defined 

by 

3%x ■= \u G ^ 2 (fi): U is constant in each interval (a, 13) C fix}- (2.18) 
Notice that we have Pji? x (V) = V a.e. in 17 \ fix and 

Pj£> x (V) = ■f V(s) ds in any maximal interval (a,/3) C fix, (2-19) 

J a 

for all V G Jz? 2 (fi). Jensen's inequality then easily yields 

Lemma 2.2 (J^x-Contraction). Let ip : R — > [0,oo] &e a convex l.s.c. function. Then for all X G we 
have 

f iP(P^ x (X))dm < / V(A)dm. 

For any pair of functions X, Y G _S? 1 (fi) we say that Y is dominated by X and we write Y -< X if the 
value of each convex integral functional on Y is less than the corresponding value on X, i.e. 



Y^A / ^(F)dm^ / ^(X)dm 



<> 



for all convex, lower semicontinuous ip: M — > [0, oo]. Estimate (12 . 16[) shows that Pj^(A) -< X for all 
AGJz* 2 (fi). 

Normal and Tangent Cones. It is immediate to check that the subdifferential (|2.1ip of the indicator function 
Jjr coincides with the normal cone of J(f at X G defined by (|1 .31|) . Applying [T71 Thm. 3.9] we 

get the following useful characterization: 

Lemma 2.3. Let X G be given. For given W G Jz? 2 (fi) we denote by 

rrn 

Ejy(m) := / W(s)ds for all m £ [0,1], 
Jo 

its primitive. Then W G Nx-^ if and only ifE>w G Afx, where 

J\fx := {s G C([0,1]): E^O in [0,1] and S = m fi\fix}- 

That is, a function 1Y is in the normal cone Nx^ if and only if it is the derivative of a nonnegative 
function S that vanishes in fi \ fix- This implies in particular that W vanishes a.e. in fi \ fix- Moreover, 
for any maximal interval (a, 0) in the open set fix we have that S(a) = S(/3) = 0, by continuity of S. Thus 

rP 

/ W(s) ds = for any maximal interval [a, p) C fix- (2.20) 

For later use, we also highlight the following fact: Let Ai, A 2 G J€ . Then 

fix, C fix 2 =>■ Ax 2 ^. (2.21) 

This follows immediately from the corresponding monotonicity for A/x . 

Let us now consider the Tangent cone TxJ^ to at A G J^: it can be defined as in (|1.18[) by 

TxJT := cl( |J 0(Jt - A)) = cl{(9(A - X) : A G JT, > o}, (2.22) 

6>>0 

or, equivalently, as the polar cone of Nx^, i.e. 

TxJT := jl7 G Jz? 2 (fi): y" ?7(m)W(m) dm < for all W € A^jr}. (2.23) 
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Lemma 2.4. Let X G Jl? be given. Then 

Tx^f = |Z7 G Jz? 2 (f2): U is nondecreasing in each interval (a, (3) C f2x;j. 

More precisely, the map U E T x J%f must be nondecreasing up to Lebesgue null sets. We may assume 
that U is right-continuous in each (a, /3) C tlx- 

Proof of Lemma \2.4\ Let U G T^X be given and fix some interval (a, /3) C 0, x . For au nonnegative 
ip G ^(fi) with spt ip C (a, /3) we have ip' G Nx-Jf because of Lemma [2731 By definition of the tangent cone 
Tx-W we find that 

rP 

/ U(m)tp' (m) dm ^ for all nonnegative </? G 3>{Q.) with spt<£> C (ct,(3). 

J a 

This shows that the distributional derivative of U in (a, /?) is a nonnegative Radon measure, and so U is 
nondecreasing in the interval. 

Conversely, assume that U G Jz? 2 (f2) is nondecreasing in each interval (a,/3) that is contained in tlx- For 
any W G Nx^ we then decompose the integral 

U(m)W(m) dm = f [/(m)T^(m) dm + V f U(m)W(m) dm, (2.24) 

where the sum is over all maximal intervals (a n , (3 n ) C £lx (at most countably many). Then the first integral 
on the right-hand side vanishes because W(m) = for a.e. m G £1 \ Qx- For each integral in the sum, an 
approximation argument (see again Lemma 3.10 in [17) ) allows us to integrate by parts to obtain 

Pn fPn ^ 

U(m)W(m)dm=- / ^ w (s)-f(ds), 



where 7 is the distributional derivative of U in (a„,/3„). Since U is assumed nondecreasing and S^/ is 
nonnegative, we conclude that U G T x Jff ■ □ 



Recalling Lemma 12.41 it is immediate to check that 

U EJ?x <=^ U and —Ue T x X '. (2.25) 

Observe that if X\ , X 2 E then 

n Xl cn X2 =^ Jf X2 CJf Xl . (2.26) 

Whenever U E Jt?x, then ([2~24|) equals zero because every term in the sum vanishes since U is constant and 
W has vanishing average. Thus 

N X X C ,^x f or all X eX, (2.27) 
with fflx the orthogonal complement of iffx ■ In particular 

Y = P Jfr (X) =^ Y = Pj^(X), ifyC*. (2.28) 

We have in fact a more precise characterization of J^x m terms of iVxJ^: in the following, let us denote 
by J*(£lx) the collection of all maximal intervals (a, f3) (the connected components) of fix- 
Lemma 2.5. For every X E X the closed subspace Jffx is 

Jfx 1 — \w E JS? 2 (f2) : W = a.e. inCl\ Q x , 

(2.29) 

W(m)dm — for every (a, (3) E ,^(Vtx) >, 

and it is the closed linear subspace of Jzf 2 (Q) generated by NxJf ■ Moreover, it admits the equivalent 
characterization 



'//' 



f = [w E^ 2 {il): I W(m)(p(X(m))dm = 0for alltpeCbiWL)}. (2.30) 
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Proof. (|2.29[) follows immediately by the definition (|2.18|) of J#x- (|2.27[) shows that the linear subspace 
generated by is contained in to prove the converse inclusion it is sufficient to check that any 

U G Jz? 2 (f2) orthogonal to all the elements of J€ is also orthogonal to i.e. it belongs to Jtf?x- This is 

true, since if U is orthogonal to Nx^ then both U and —U belongs to the polar cone to Nx^ which is 
T x J(f: by ([235)1 we deduce that U G 3%x- 

Concerning (|2.30|) we simply notice that all U G J4?x can be written as U = uoX for a map u G J?? 2 (R, g), 
where g = X#m. Approximating u in L 2 (R, g) by a family of functions <p G C(,(M) we obtain (12.301) . □ 

Lemma 2.6. For any X G and U G Tx^ we have that 

(P^ x - id)U G N x .Je. (2.31) 

Proof. Lemma 3.11 in [J7] shows that (TOT]) holds if U G J^. Since Pj*> x X - X = 0, (|2~BT1) holds for f7 - X 
and, since Nx<$f is a cone, for arbitrary 0(t7 — X), 9 > and [7 G J^. We conclude recalling (|2.22l) . □ 

Remark 2.7. If ^ is a closed convex subset of if 2 (ft) and X G L 1 ((0,T);X) with G ^ for a.e. i G (0,T) 
then it is easy to check that 

-f X(t)olt G ^ ; moreover f X{t)dt^ c g if <^ is a cone. (2.32) 
Jo Jo 
In fact, applying Jensen's inequality to the indicator function of ^ we get 



X{t)&b)<j I v (X(t))dt = 



If in addition c <§ is a cone then I<%(TX) = I<g(X) for every X G Jz? 2 (f2) and we deduce the second implication 
of (I2~32l) . 

2.3. Convex Functions. In this section we recall some auxiliary results on convex functions. We are 
interested in functions ip : R — > [0, oo) that are 

even, convex, of class C 1 (R), with ip(0) = 0, (2.33) 

and for which the homogeneous doubling condition holds: 

there exists q > 1 such that ip(\r) < \ q ip(r) for all r G R, A ^ 1. (2.34) 

Notice that if condition (|2.34l) holds for ip, then it also holds for the map r H> ip 2 (r), with exponent qp. 
Combining (|2.33p and (I2.34[) . we obtain the inequality 

ip(ri +r 2 ) sC 2 q - 1 (ij;(r 1 ) + ip{r 2 )) for all n , r 2 G R. (2.35) 
We will denote by ^ : «5f 1 (J7) — > [0, oo] the associated convex functional 

:= / ?/>(X(m))dm for all X G Jzf^fi). (2.36) 
Jn 

Lemma 2.8. Suppose ip: R — > [0, oo) satisfies (|2.33[) . TTien </ie doubling condition (|2.34[) Zio/ds if and only 
if ip has one of the following, equivalent properties: 

rip'(r) < #(r) /or a^Z r > 0; (2.37) 

i/iere exists C ^ suc/i i/iai 0(2r) < CV>(r) /or aZZ r > 0. (2.38) 

Proof. Property (|2.38p is a consequence of (|2.34l) . and (12.371) follows from 

,u \ ,. i>(\r)-ip(r) A 13 — 1 . . 

r^> (r) = hm $C hm — VA r ) — QVv) 

a^i+ A — 1 A-S-1+ A — 1 

for all r > and q > 1. 

To prove the converse statement, we notice first that since "0 is an even, smooth function, we have that 
ip'(0) = and so ip is nonnegative and nondecreasing for all r > 0, by convexity. Moreover, if (|2.38[) holds, 
then again by convexity we find 

W< ttM <(C _i)W forallr>0. 
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Thus (|2.37l) holds with q :— C — 1, which must not only be a nonnegative number but must be greater than 
or equal to 1. 

Assuming now that (|2 .37[) is true, we consider the Cauchy problem 

i]'(s) — q for s G [r, oo), with r?(r) = ip(r), (2.39) 

s 

which admits a unique solution 77(5) = ip(r){r~ l s) q for all s r > 0. A standard comparison estimate for 
solutions of ordinary differential equation yields 

ip(s) < r?(s) = (-)VO - ) for all s ^ r > 0. (2.40) 

Since ip is nondecreasing, we conclude that q ^ 1 and then ()2.34p follows for r > 0. By evenness of ip and 
since ip(0) = 0, the inequality extends to r ^ as well. □ 

Lemma 2.9. Let p G [1, 00) 6e given and supppose that the function n : R — ^ [0, cx)) satisfies (|2.33[) and </ie 
p-coercivity condition 

0<liminf^V lim ^ = 00. (2.41) 

r— >0+ r r— >oo r 
For every q > p, there exists a map ip: R — > [0, 00) satisfying (|2.33[) / (I2.34p wii/i 

?/>(r) < ?j(r) /or aZI r 6 M and lim = 00. (2.42) 

Proof. By |191 Lemma 3.7] it is not restrictive to assume that 77 is of the form fj 2 for a suitable convex function 
fj with superlinear growth, and so we may just consider the case p = 1 (see the remark following (|2.34p ). By 
convolution, we can assume that n is smooth in the open interval (0, 00), with 5 := inf r> o r]'(r) > 0. 
We then choose g > 1 and we set ip( r ) '■— Sr q /q for all r G [0, 1], so that 

tp(r) < r?(r) in [0, 1] and ^'(1) = 8 < r/(l). 

For r^lwe define ^ to be the solution of the Cauchy problem 

tp'(r) =mhi\r]'(r),q^-^ 1 for r G [1, 00), with = S/q. (2.43) 



Then ip(r) ^ 77(f) for all r ^ and ?/> satisfies (|2.37p of the previous lemma. 

To prove that ip also satisfies (|2.33[) . notice first that ip'(0) = since q > 1, and that ip' is continuous at 
7- = 1. Hence ip can be extended to an even C 1 (M)-function. In order to check that ip' is nondecreasing, let us 
first observe as a general fact that if a continuous function /? is nondecreasing in each connected component 
of an open set icl that is dense in [1, +00), then [3 is nondecreasing in [1, +00). We apply this observation 
to (3 :— ip' and we set A := Aq U A\, where 

A, :=/re(l,oo):?^<»/(r)j, 



A\ := interior of ^ r G (1, 00) : 77' (r) ^ q 



In each connected component of Aq, the function ip solves the differential equation ip'(r) = qip{r)/r, and so 
ip is of the form cr q for a suitable constant c > 0. Therefore ip' is nondecreasing in A . On the other hand, 
in each connected component of A\, we have that ip'{r) — rf(r) and r( is nondecreasing, by assumption. 
Finally, notice that ip is nondecreasing on the interval [0,1] since ip(r) = 5r q /q there. We can now apply 
Lemma [^U to conclude that ip has the doubling property (I2.34j) . 

It only remains to prove the second statement in (|2.42[) . Since 77 has superlinear growth, its derivative 
77' (r) — > 00 as r — >• 00. Assume now that ip(r)/r remains bounded as r — > 00. Then there exists a number 
7*1 ^ 1 such that 

ip'{r) = q for all r G [ri,oo); 

r 

see (|2.43|) . But this implies that ip'(r) — cr q + co for all r G [ri, 00) and suitable constant c > and therefore 
is unbounded as r — > 00. This is a contradiction. □ 
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Lemma 2.10 (Compactness in J(f). Let ^ be the integral junctional defined in (|2.36[) corresponding to an 
even, convex function ip: R — > [0,oo) with 

lim Pw = oo. (2.44) 
|rK«, |rf 

Then each sublsevel of \1/ 

X(^,a) := {X G X: $[X] < a} is compact in ^f 2 (fl) for every a>0. 

Proof. Because of (|2.44|1 . the _Sf 2 (O)-norm of elements of J(f(^,a) is bounded by some constant A that 
depends on a and ^ only. By monotonicity, we find that 

1 r 1 fir 1 x 1/2 

X(w) < / A(m)dm< / |A(m)| 2 dm 

(1 — w) 1 /^ 

for all X G a). Analogously, we obtain a lower bound 

X(u>) ^ — - for all w G Vl. 

Any sequence {A n } in JT(\E', a) is therefore uniformly bounded in each compact interval [8, 1 — S] where 
5 > 0. Applying Helly's theorem and a standard diagonal argument we can find a subsequence (still denoted 
by {X n } for simplicity) that converges pointwise to an element X G JXf. Since ip satisfies ()2.44j) . the sequence 
{\X n — A| 2 } is uniformly integrable and thus X n — !• X in jSf 2 (O) . □ 

3. Lagrangian solutions 

As explained in the Introduction, when studying system (jl.lj) . one is lead to consider solutions to the 
Cauchy problem for the first-order differential inclusion in Jz? 2 (f2) 

X{t) + dI,ytf{X(t)) 3 V + f F[X(s)] ds for t ^ 0, X(0) = lim X(t) = X, (3.1) 

Jo *4-0 

and, possibly, satisfying further properties. 

Before discussing (|3.ip . we will state below the precise assumptions on the force operator F; examples, 
covering the case of (|1.2p or (|1.3[> . are detailed in Section [51 

3.1. The force operator F. Let us first recall the link of the map F: Jif — >• „S? 2 (f2) with the force 
distribution /: £P(R) — > J( (R) in (pTTj) : as in ([04"]) we will assume that 

ip(x) f[g](dx) = I tp(X e (m))F[X](m)dm for all tp G g G &>(R) (3.2) 

Jn 

recalling (|2.30[) one immediately sees that F[X] is uniquely characterized by (|3.2I) only when = {0} or, 
equivalently, when J%x — Jz? 2 (fl) i.e. fix = 0: this is precisely the case when X is (essentially) strictly 
increasing. 

One could, of course, always take the orthogonal projection of F[X] onto J#x in order to characterize it 
starting from (I3.2[) . This procedure, however, could lead to a discontinuous operator which would be hard to 
treat by the theory of first order differential inclusions. This happens, e.g., for the (attractive or repulsive) 
Euler-Poisson system. We thus prefer to allow for a greater flexibility in the choice of F complying with (|3.2j) , 
asking that it is everywhere defined on Jf and satisfies suitable boundedness and continuity properties. 

Definition 3.1 (Boundedness). An operator F: X — > Jz? 2 (f2) is bounded if there exists a constant C ^ 
such that 

||F[X]||^ (n) <C 2 (l + pf|U« ( n)) for all AG JT. (3.3) 
We say that F is pointwise linearly bounded if there exists a constant C p ^ such that 

\F[X}(m)\ ^C p (l + \X(m)\ + \\X\\^i m ^ for a.e. m G fl and all X G X. (3.4) 
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Note that if F is pointwise linearly bounded, then F is bounded and satisfies (I3.3[) with the constant 
C 2 := 2C p . 

Let us recall that a modulus of continuity is a concave continuous function u>: [0, oo) — > [0, oo) with the 
property that = u>(0) < u(r) for all r > 0. 

Definition 3.2 (Uniform continuity). We say that an operator F: X — > Jzf 2 (£l) is uniformly continuous 
if it is bounded as in Definition 13.11 and there exists a modulus of continuity uj with the property that 

llF^-i^H^^ <w(||Xi-* 2 |Ur*(n)) for all X 1; X 2 e JT. (3.5) 

We say that F is Lipschitz continuous if it is uniformly continuous and (|3.5j) holds with u){r) — Lr for all 
r 0, where L is some constant. 

Notice that if F is uniformly continuous then it is also bounded. Whenever a uniformly continuous F is 
defined by (|3.2|) on the convex subset J^i of all the strictly increasing maps and satisfies (|3.5j) in J€ai, then 
it admits a unique extension to J?T preserving the continuity property (|3 . 5|) and the compatibility condition 

(EUD . 

As we observed at the beginning of this section, a last property of F which will play a crucial role concerns 
its behaviour on the subset Qx where the map X is constant. Since the force functional determines the 
change in velocity, in the framework of sticky evolution it would be natural to assume that 

F[X] G Jfx for every XeX. 

We shall see that a weaker peroperty is still sufficient to preserve the sticky condition: it will turn particularly 
useful when the attractive Euler-Poisson equation will be considered. 

Definition 3.3 (Sticking). The map F: X — > Jz? 2 (£l) is called sticking if for all transport maps I,i"e Jff 
with Y G J4?x we have 

F[X]-P^ x (F[X])edI^(Y). 

3.2. Lagrangian Solutions. Let us start by giving a suitable notion of solutions to (|3.1[) . 

Definition 3.4 (Lagrangian solutions to the differential inclusion (I3.1[) L Let F : -J€ — > L 2 (Vl) be a 
uniformly continuous operator and let X G 3€ and V G 3%? — jSf 2 (O) be given. A Lagrangian solution 
to (|3.ip with initial data (X, V) is a curve X G Lip loc ([0, oo); J(f) satisfying X(0) = X and (|3.ip for a.e. 
tG(0,oo). 

By introducing the new variable 

Y(t) :=V+ I F[X(s)]ds 
Jo 

we immediately see that (14. ip is equivalent to the evolution system 

f X(t) + dIjr(X(t)) 3 Y(t), _ _ 

^ fort>0, 10,70 = (X,V), (3.6) 

\ Y(t) = F[X(t)}, 

Notice that the continuity of F yields Y G C 1 ([0, oo); _Sf 2 (SI)) . We state in the following Theorem the main 
properties of the solution X to (|3.ip 

Theorem 3.5. Let F : „S? 2 (f2) — > Jtf be a uniformly continuous operator and let (X,Y) be a solution to 
(|3.6p . Then the following properties hold: 

• Right-Derivative: 

d+ 

The right-derivative V := "^~A exists for all t ^ 0. (3-7) 

• Minimal Selection: 

V(t) = (Y(t) - dIjr(X(t))Y for all t > 0, (3.8) 

where A° denotes the unique element of minimal norm in any closed convex set of A C Jz? 2 (f2). In 
particular if we replace X(t) by V(t) then (|3.1I) and (|3.6p hold for all t ^ 0. 

• Projection on the tangent cone: 

V{t) = P Txm *r{Y{t)) far ail t (3.9) 
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• Continuity of the velocity: 

V is right- continuous for all t ^ 0; (3.10) 

in particular 

limF(t) = V if and only if V 6 T x JfT. (3.11) 

If T a C (0, oo) is the subset of all times at which the map s \— > \\V(s)\\^2^qj is continuous, then 
(0,oo) \ T° is negligible and at every point of T° V is continuous and X is differ entiable in _Sf 2 (f2). 
Setting g t := X(t)#m there exists a unique map Vt £ £ 2 (R, gt) such that 

X(t) = V(t) = P^ x(t) (Y(i)) =v t oX t e Jf x(t) for every t 6 T° . (3.12) 

• Solution to (|l.lj) : If moreover F is linked to f by (|3.2j) . g = X#m and V — v o X, then the couple 
{g, v) defined as above is a distributional solution to such that 

]img(t,-) = g in ^ 2 (K), lim g(t, -)v(t, •) = gv in JK(R). (3.13) 

Proof. p.7|) . p.8|) . and f|3.10|) are consequence of the general theory of [8], Theorem 3.5; ()3.9|) follows 
immediately by §^ since V(t) 6 T x(4 )Jf by (JSlTj) and V(t) + N X (ty^ 3 ^(i)- 

Concerning (|3.12l) we can apply the Remark 3.9 (but see also Remark 3.4) of [5], which shows that at 
each differentiability point i of X its derivative is the projection of onto the affine space generated by 
Y(t) — dIje(X(t)), i.e. the orthogonal projection of Y(t) onto the orthogonal complement of the space 
generated by dI^(X(t)). Recalling Lemma |2~B1 we get (|3.12l) . 

In order to prove the last statement, we use the crucial information of (|3.12l) that V(t) £ J#x(t) f° r a - e - 
t > 0, a fact which may have been noticed for the first time in|llj. In particular, the projected velocities 

V{t)=P jexm V{t) t>0, (3.14) 

concide with V(t) for every t G T°, where T° is a set of full measure in (0,oo). Since any element V of 
J^x(t) can be written as v o X for a suitable Borel map v e J^ 2 (il) we deduce that there exists a Borel map 
v : [0, oo) x M ->• R such that v(t, •) £ ^f 2 (M, g(t, ■)) and 

V*(t,-) =v(t,X{t,-)) a.e. in O, for every t^s 0. (3.15) 

From Equation (13. 14)) we also have V(t) = v(t,X(t)) for t € T°. We then argue as follows: For all test 
functions ip £ ^((0,T) x R) we have 



dtsp(t, x)v(t, x) + d x ip(t, x)v 2 (t, x)j g(t, dx) dt 

d t ip(t, X(t, rn)) + d x (p(t, X(t, m))v(t, X(t, m)) ) V(t, to) dm dt 



d t tp(t, X(t, to)) + d x ip(t, X(t, m))v(t, X(t, m))j Y(t, to) dm dt 

— Vpit, X(t, m))l ) Y(t, m) dm dt 
dt J 

<p(t, X(t, m))F[X(t, -)](to) dm dt (3.16) 
Applying formula (|3.2[) in ()3.16[) . we obtain 

dttp(t, x)v(t, x) + d x (p(t, x)v 2 (t, x) I g(t, dx) dt 



ip(t,x)f[g(t,-)}(dx)dt, (3.17) 

which yields the momentum equation in (|l.lj) in distributional sense. An (even easier) analogous argument 
holds for the continuity equation. This shows that the pair (g, v) defined by (|3.14[) and (|3.15[) is a solution 
of (TTI]). 
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The first limit of (|3.13l) follows since \im t ^ X(t) = X in J2? 2 (S1) and X#m =: g. Concerning the second 
limit of (|3 . 13[) we have to show that 

<pv(x)g(dx) = lim / ipv(t 7 x)g(t,dx) (3.18) 



£2 



for every ip £ C 6 (R). Since V = v o X we have V £ M" x C 7>Jf so that lim t;o V(t) = V in «5f 2 (ft) and 
therefore 



V?(x)v(x)p(dx) = / (p(X)Vdm = Km / <p(X(t))V(t) dm 



= lim / <p(X(t))V*(t) dm = lim / (p(x)v(t,x)g(t,dx), 
where we used the fact that V{t) — V*(t) is perpendicular to Jfx(t)- D 
As we already observed in the previous proof, notice that (|3.11[) surely holds if V € . 

3.3. Existence, uniqueness, and stability of Lagrangian solutions for Lipschitz forces. Applying 
the general results of [3] is not difficult to prove 

Theorem 3.6. Let us suppose that F : JT -> Jz? 2 (ft) is Lipschitz. Then for every (X,V) e x Jz? 2 (ft) 
i/iere exists a unique Lagrangian solution X to (|3.1[) and /or every T £/iere exists a constant Ct 
independent of the initial data such that for every t € [0, T] 

||*(*)||j^(Q) + ||V(*)||js?'(n) < C T (l + ||A||^ (0 ) + ||V|U«(n)). (3.19) 

Moreover, for any T ^ i/iere exists a constant Ct wi/i £/ie following property: For any pair of strong 
Lagrangian solutions Xi with initial data (Xi, Vi) for i = 1, 2 we /lave i/iai /or t £ [0,T] if Zio'ds 

||Xi(t) -X 2 (t)||^ 2(n) < Ct(\\Xi ~ X 2 \\sf*w + II VI - V 3 \\s»(a)) (3-20) 
/ || Vi(t)-V 2 (t) \\% 2{n) dt^C T J2 (ll**ll*»(n) + IMUf'(n)) 

■'° ..2 

x (||A > 1 -A > 2 ||^ 2(n) + -Kalian)). (3.21) 

Proof. Recalling the equivalent formulation (|3.6|) . we introduce the Hilbert space H : Jz? 2 (f2) x L 2 (fl) and 
the (multivalued) operator A(X,Y) :— (dlj^(X) — Y,F[X]). It is easy to check that A is a Lipschitz 
perturbation of the subdifferential of the proper, convex, and l.s.c. functional Y) := I^(X). Thus 
existence, uniqueness, and the estimates (|3.19p . (I3.20[) follow by [HI Theorem 3.17]. 
The same estimate also yields for the second component Yi(t) = V + L F[Xi(s)] ds 

\\Yi(t) - Y 2 (t)\\^ {n) < C T (\\X X - X 2 \\se* m + \\V - V 2 \\^ {n) ) (3.22) 
and, by the boundedness of F, 

f \\Yi(t)\\%*(n)to < c T (i + \\x t \\% 2{n) + 11^11^,(0)). 

Applying Theorem 2 in [20 to the first equation of (|3.6[) we get (|3.21[) . □ 



A straightforward application of the previous Theorem shows that Lagrangian solutions are stable if F 
is Lipschitz: a sequence of Lagrangian solutions with strongly converging initial data converges to another 
Lagrangian solution. 
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3.4. Sticky lagrangian solutions and the semigroup property. We consider here an important class 
of Lagrangian solutions. 

Definition 3.7 (Sticky lagrangian solutions). We say that a Lagrangian solution X is sticky if 

for any t\ ^ t 2 we have Q,x(ti) C Clx(t 2 )- (3.23) 
By (|2.2ip and (|2.26p any sticky Lagrangian solution satisfies the monotonicity condition 

0/jr(*(ti)) C dIjr(X(t a )), C J^x( tl ) faranyti<i 2 . (3.24) 
The nice features of sticky Lagrangian solutions are summarized in the next results. 
Proposition 3.8 (Projection formula). If X is a sticky Lagrangian solution then 

V(t) G Mx(t) for all times t ^ (3.25) 

and it satisfies 

X(t) = Pjr (x + j Y{s) ds\ = Pj<r (x + tV + J (t - s)F[X(s)] ds^j , (3.26) 

V(t) = Pje xw (Y(t)) = P^ xm (? + j F[X{s)] ds) . (3.27) 

Proof. (|3.25|) follows easily by (|3.10p and (|3.12p , thanks to the monotonicity property (|3.24p . Equation (|3.27l) 
then follows from Q3.6P (where X is replaced by V) and (|2.27D . 
In order to prove (|3.26p . we set for s ^ 

E(s) := Y(s) - V(s) G dI. x {X{s)) C dI x {X(t)) Ht>s, (3.28) 

and we integrate (|3.28p w.r.t. s from to t to obtain 



E(s)ds = X - X(t) + / Y(s)ds G dIjr(X(t)) 
Jo 

Recalling (l2~T3|) we get (j3~26|) . □ 

Lemma 3.9 (Concatenation property). Let X\,X 2 be Lagrangian solutions with initial data X\,V\ and 
X2, V 2 respectively and let us suppose that 

Q X2 C n X2(t) for every t > 0. (3.29) 

If for some t > 

X 2 = X 1 (t), Y 1 (r)-V 2 =E 2 edl r (X 2 ), (3.30) 

then the curve 

X:=f X f s 1/0 (3.31) 

\X 2 (t-T) ift^T, V ^ 

is a Lagrangian solution with initial data (X\,Vi). In particular, if X\,X 2 are sticky Lagrangian solutions, 
then X is also sticky. 

Notice that 

the choice V 2 := V{t\) always satisfies (|3. 301) . (3.32) 
Proof. It is easy to check that 

Y (t)-=i Yl{t) ifO^t^r, 
U ' \Y 1 {T)-V 2 +Y 2 (t-T) ift>r, 

is Lipschitz continuous and satisfies 4fY(t) = F[X(t)] a.e. in (0, 00). We have to check that X satisfies the 



first differential inclusion of (|3.6p for t ^ r w.r.t. Y. By definition of X we have for t ^ r 

Y(t) - ix(t) = Y 1 (t) ~V 2 + Y 2 (t - t) - V 2 (t - t) 
dt 

= E 2 + E 2 (t -r) G dI x (X 2 {t - r)) = dIjt(X(t)) 
since by (|3"3tjj) and (j3"^9l E 2 G dIjr(X 2 (t - r)). □ 
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It would not be difficult to show that Lagrangian solutions in general do not satisfy the sticky property 
nor the semigroup property. If the force is sticking then the next property shows that these properties are 
strictly related. 

Theorem 3.10 (Semigroup property). If the force operator F is Lipschitz and sticking and 

Q x C fijf(t) f or every Lagrangian solution X starting from (X, V) G ,X x ,¥P X (3.33) 

then every Lagrangian solution X starting from (X, V) G Jff x ,¥tf x is sticky and satisfies the following 
semigroup property: for every r > the curve X{t) := X(t — t) is the unique Lagrangian solution with 
initial data X(t), V(t). 

In particular, for all t ^ t\ > we have 

X(t) = ?. x (x(t x ) + (i - tjVfa) + (t - s)F[X{s)} ds^j (3.34) 

V(t) = P^ xw (v(h) + J* F[X(s)] ds) . (3.35) 

Proof. Let T° G [0,oo) as in (13.121) ((0,oo) \ T° is negligible). For every r G T° consider the Lagrangian 
solution X2 with initial datum (X(t), V(t)): by the concatenation property (with the choice p. 321) ) the 
map X defined as in p. 311) (with X\ := X) is a Lagrangian solution and therefore coincides with X, since 
F is Lipschitz. (|3.33[) yields that 

Hx(t) C n x{t) for every O^t <t, t G T° U {0}. (3.36) 

Let us now fix s > and consider a sequence h n J, such that 

-L(X(s)-X(s-h n )) mL 2 (n). 

fori 

Since T x ^J(f is a closed convex cone, it is also weakly closed, so that by its very definition definition we 
have -V- G T x{s) X . 

We set S(i) :— Y(t) — X(t) G dIj^(X(t)) thanks to the differential inclusion of (|3.6p ; an integration in 
time from s — h n to s and (|3.36[) yield 

Y (r) dr - ^-(X(s) - X(s - h n )) = / 3(r) dr G dI, x {X{s)). 

Passing to the limit as n — ¥ 00 we obtain 

£_ :=Y(s)-V-e8Ijr(X(8)) 

and therefore by (|2.27|) 

V:=Px xw (Y-) = P XxM (Y(s)). 

Since 

Y( S ) -v = y-V- + (-v- (-v-)) = »_ + (-?- (-y_)) g a/jrWs)) 

by (I2.31[) and the fact that — V- G T x ^ s y¥f , we can apply the concatenation property as before, joining at 
the time s the Lagrangian solution X\ := X with the Lagrangian solution X 2 arising from the initial data 
X := X(s) and V. The uniqueness theorem shows that this map coincides with X and therefore (|3.29[) 
yields ttx( s ) C fijt(t) f° r every t > s. 

In particular we have V(t) G ^x(t) for every t > so that a further application of the concatenation 
Lemma 13.91 yields the semigroup property. (I3.34[) and (|3.35[) follow then by the corresponding (|3.26p and 

(E2ZD ' □ 

We conclude this section with our main result conerning the existence of sticky Lagrangian solution; the 
proof will require a careful analysis of the discrete particle models and therefore will be postponed at the 
end of Section [SJ see Remark 15.41 

Theorem 3.11 (Sticking forces yields sticky Lagrangian solutions). // the force operator F is Lipschitz and 
sticking (according to Defintion \3.3\) then every Lagrangian solution to (|3.1[) with X G and V G Jt? x is 
sticky. 
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Remark 3.12. We have seen that the right-derivative V of a sticky Lagrangian solution is right-continuous 
everywhere. It is continuous for all t ^ for which the function t \-> \\V(t)\\^2^Q-j (which represents the 
kinetic energy) is continuous; see Proposition 3.3 in [5]. At such times the map t h->- X(t) is differentiable. 
We do not know whether the velocity is of bounded variation. But (|3.35p and (13.91) show the following 
statement: For any t ^ let V_ £ Jz? 2 (f2) be any weak accumulation point of V(s) as s j" t. Then 
V(t) — P Tx jg-{V_) = P^ x(t) (V—). This is the analogue of the impact law (|1.20|) we discussed in the 
Introduction. It follows easily from (13 . 35[) . 

3.5. Lagrangian solutions for continuous force fields. The goal of this section is to extend the existence 
Theorem 13.61 to the case of (uniformly) continuous force operators. 

Theorem 3.13. Suppose that F : J£ — > J^f 2 (Q) satisfies the pointwise linear condition \3.$ and it is 
uniformly continuous according to V3. 5\) . Then for every (X, V) £ ,Jtif x Jzf 2 (Q) there exists a Lagrangian 
solution (X,Y) of |[3T5jl . 

Moreover, for any T there exists a constant Ct ^ such that any Lagrangian solution X with velocity 
V := ^X satisfy 

ll^f 2 (n) + ll^ll^ 2 (o)) (3.37) 

for all t £ [0,T]. Iftp-.M. — > [0, +00) is an integrand satisfying (|2.33|) and (|2.34p for some q ^ 1, then there 
exists a constant C q ,r such that 

9[X(t)] + *[V(t)] 5$ C q:T (l + V[X] + *fV]) (3.38) 

for all t £ [0, T], with functional ^ defined in (|2.36p . 

Proof. It suffices to show that there exists a solution to (|3.6[) in a bounded interval [0, T) with T independent 
of the initial condition. We will choose 

1 

T := < 1. 

2C P 

where C p is the constant of (|3.4p . which is not restrictive to assume greater than 1. 

We consider the following operators defined in C([0, T]; «5? 2 (£1)): the first one O x maps W £ C([0, T]; Jz? 2 (£l)) 
into Y := 0\{W) defined by 

Y(t):=V+ [ F[W{s)]ds; (3.39) 



the second one, 2 , maps Y £ C([0, T]; L 2 (Vl)) into the solution X — 2 {Y) of the differential inclusion 

X + dL. x (X)3Y, X(0)=X. (3.40) 

Both of them are continuous, since 

||0i(Wi) - 2 (W 2 )\\ 0C < Tw(]|Wi - Walloo) (3.41) 

(where we denoted by || ■ ||oo the usual sup norm in C([0, T]; Jzf 2 (f2))) and 

\\O2p1) - 2 (Y 2 )\\ 00 < T||Y! - Y 2 \\ oa . (3.42) 

We want to show that O := 2 o 0\ has a fixed point X, which is a Lagrangian solution with initial data 
(X,V). We may use de la Vallee Poussin Theorem and Lemma [2.91 to obtain ip: R — > [0, 00) satisfying 
P35|l / ([23SP for some q > 2 and (using the notation ([235]) ) 

Um ^1 = 00, [X] + < 00. (3.43) 



Choose m large enough so that 



(11 



and let 
where 



|y|| 2 + 2C p (l + m+||y|| 2 )JT<m, (3.44) 
JT(tf) = {VK e JT: ||W-X||^.2 (n) < m, v(w-X) < D} 
A(T) + (12CpT)«T^(m) + 
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and 

A(i) := (t9(V) + (3T)«- 1 C 9 ^(l)y ) . 

We eventually set 

^ = {l£ C([0,T];J^ 2 (ft)): X(t) e JXT(^) Vt G [0,T], 

||X( S )-X(r)||^ ( o) < -|r- S | Vr.se [0,T]j, 

which by Arzela-Ascoli Theorem is a nonempty, compact, and convex subset of C([0, T]; Jz? 2 (i!)). In light of 
the Schauder Fixed Point Theorem it suffices to show that O maps 'io into itself. 

Let W € and set Y = O x {W),X = 2 (Y) = 0{W) with V = ^X. We exploit Lemma EH and 
equation (|3.53|) below with Jensen's inequality to obtain 



V(x(t)-X) = V(J V(s)ds 



< T^ 1 

< T^ 1 



A(t) + (6C p T 2 y / *(W(l))dl 



A(t) + (12C p r)9 y (*(W(i) - X) + 
A(t) + (12C p T)T(V(m) + 



< D. 



(3.45) 



Using ^(r) = r 2 in (|3 .48[) we have 



\\V(t)h < \\V+ / F[W(s)]dfl|L < ||V-|| a + / ||F[W(s)]|Lds 



where we have used the J^x -Contraction property in Lemma 12.21 We use that F is also bounded with 
constant C 2 = 2C p and that Y S ^ to conclude that 



||F(t)|| 2 <||y|| 2 + 2C p (l + |m| + ||A|| 2 )<-. 



(3.46) 



Thus, for every ^ r ^ s ^ T 



||X(r) - X( S )|| 2 < J \\V{s)\\ ds < -\r - s\, \\X(s) - X\\ 2 < m. 

These prove that 

Concerning the estimates (|3.37[) and (|3.38p we simply set W = X in the next Lemma 13.141 and apply 
Gronwall's lemma. □ 

We conclude this section with the uniform bounds for solutions to differential inclusions invoked by the 
previous fixed point argument. In view of the next applications, we state them in a slightly more general 
form. 

Lemma 3.14 (A priori bounds). Let F: J(f — > „S? 2 (f2) be pointwise linearly bounded so that there exists 
C p > such that Jff.^| ) holds. Let ijj : M. — 5- [0, +oo) be an integrand satisfying (|2.33[) and (|2.34p for some 
q^l. LetX € Lip(0, T; ^ 2 (0)) and Z,W <E JC°°(0, T; Jz? 2 (£l)) be such that X (0) = X and 

V(t)+dIjt(X(t)) 3 V+ [ Z(s)ds Vte[0,T) (3.47) 



where 



Then for a.e. t £ (0, T) 



V{t) 



d+X 

dt 



(t) and Z(s) ~< F[W(s)]. 



#(V(t)) < * [V+ / Z(s)ds 



(V(t)) < 2?- 1 U(V) + {itf- 1 ^ (i>{l)t + 2 f *(W(a))<L 



(3.48) 
(3.49) 
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and 



#(X(t)) < 2 9 " 1 U>{X) + T q A(T) + (6C P T 2 )« j ^{W(s))ds^j (3.50) 



Proof. Recalling Theorem 13.51 and (|3. 121) . equation (|3.47l) yields 



V(t) = Pje xlt) (V + f o Z{s)ds^j for every teT°, 



where T° has full measure in (0, T). Hence, by Lemma [2~2l 

ft , \ , rt 



which proves (|3.48|) . 

Using (|2 .35|) and Jensen's inequality we obtain 



*(V(t)) < (*(U) + 9(f Z(s)ds)^ < V- 1 (%(V) + ^ *(Z(«))ds) . (3.51) 



We use the fact that F is linearly bounded, ip is even, and VKII^IIi) ^ $(W) by Jensen's inequality, to find 
to obtain that for all s ^ 

< ¥[F(W(»))] < tf[C p (l + + H^llx)] < 3«- 1 C«(^(l) + 2*[W(*)]). (3.52) 

The first inequality in Equation (|3.52[) was obtained via Lemma 12.21 We combine Equations (|3.51[ 13.52[) to 
obtain Equation (jSMl . By Equation (1X15)) 



L 



f V(V(s))ds < 29- 1 ^(v) + (3r)*- 1 q(?A(i)y + 2^ J ^{W{l))di^j 



= A(t) + 2(6T) ? - 1 C p 9 / (i - /)*(W(0)«K 







where 



We have 



< A(i) + (6TC P ) 9 / *(W(0)tB (3.53) 
A(i) := 2'- 1 (W(V) + (3T)«- 1 C7 p ^(l)^) 



*(Jf(t)) = ^[X + V{s)ds^ <2 q - 1 [^{X)+t q - 1 V(V{s))dsj 
where have used (|2.35l) and then Jensen's inequality. This, together with (|3.53[) yields (13.501) . □ 

4. The semigroup property and generalized Lagrangian solutions 

We have seen that Lagrangian solution may fail to satisfy the semigroup property in the natural phase 
space for the variables (X, V), V = X (stated in Proposition p.lOj) for sticky Lagrangian solutions). In fact, 
the formulation given by the system (|3.6[) shows that the natural variables for the semigroup property are 
the couple (X, Y). 

This motivates an alternate notion of solution (still linked to (jl.lj) ) which tries to recover a mild semigroup 
property, at the price of loosing uniqueness with respect to initial data. 

Recall that for any transport map X € J€ the orthogonal projection Pj^ x onto the closed subspace 
C Jz? 2 (£!) leaves the given function unchanged in ft \ ilx and replaces it with its average in every 
maximal interval (a,j3) C Qx', see (|2.19p . As a consequence, the function P,^ X (F[X]) is constant wherever 
X is. 



Definition 4.1. A generalized solution to (|3.1I) is a curve X G Lip loc ([0, 00); ,Jtf) such that 
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(1) Differential inclusion: 

X(t) + dIje{X{t)) 3 V + [ Z(s)ds for a.e. t G (0, oo), (4.1) 

Jo 

for some map Z € L^ c ([0, oo); .5f 2 (ft)) with 

Z-F[X(t)]eJf?x(t) and Z-<F[X(t)] for a.e. t 6 (0, oo). (4.2) 

(2) Semigroup property: For alH ^ ii ^ the right derivative V = ^jr^ satisfies 

V(t) + dIjt(X(t)) 3 V(h) + f Z{s) ds, (4.3) 

(3) Projection formula: For all t ^ ti ^ 

X(i 2 ) = ^(tx) + (t 2 - ti)V(ti) + ^ \t 2 - s)Z(s) ds) , (4.4) 

Note that for generalized Lagrangian solutions the semigroup property and the projection one (|4.4p 
are part of the definition, while for sticky Lagrangian solutions (|3.34[) and (|3.35|) are consequences of the 
monotonicity property (|3.23[) . The obvious choice in (|4.2I) is Z(t) := F[X(t)\ for all times t ^ 0, which also 
shows that any sticky Lagrangian solution is a weak solution. 

Remark 4.2. If one is ultimately interested only in the existence of solutions to the conservation law for 
this purpose any Z stisfying (|4.2I) is sufficient. In fact, we proved in Theorem 13 . 51 that if the force functional 
F[X] is induced by an Eulerian force field f[g], so that (|3.2[) holds whenever X £ J?T and X#(m) = g, then 
any strong Lagrangian solution yields a solution of the conservation law The same argument works 

for weak Lagrangian solutions. Because of (|4.2I) we have that Pj#> x(t) (Z(t)) = Pj#> x(t) (F[X(t)]). On the other 
hand, it holds 

r ip(X(m))F[X](m)dm = f ip{X(m))P,^ x (F[X])(m) dm 



for all ip £ f^(R), with a similar formula for Z in place of F[X\. Then the argument on page [20] can be 
adapted to prove the claim; see in particular (|3.17p . 

Since X is everywhere right differentiable, we have V(t) S T X (t)^ f° r every t 0, so that (|4.3p yields 
V(t) = P Txm Jtr (v(h) + jf 2 FLY(s)] ds) for every t> h > 0, (4.5) 

which also yields 

V(t) = Pje xw (v(ti) + £ F[X(s)}ds^j for almost every t ^ ti > 0. (4.6) 

It is immediate to check that any solution is also a generalized solution, corresponding to the choice H(t) := 0. 
By introducing the new variable 

Y(t) :=V+ [ Z(s)ds 
Jo 

we easily see that (|4.1I) is equivalent to the evolution system 

(X(t) + dIjr{X{t))3Y(t), 

Y(t) = Z(t), - - 

Y' fart>0, =1,7 (4.7) 

Z(t) -< F[X(t)] 
where i? = in the case of (|3.1[) . 
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4.1. Stability of generalized Lagrangian solutions. In this section, we will prove a stability result 
for generalized Lagrangian solutions. Instead of relying on a semigroup estimate, strong compactness now 
follows from an argument based on Helly's theorem (recall Lemma fe.lOp and on the closure properties of the 
map X h-> Pjr x (F[X]) for lei. 

Lemma 4.3. Consider {(X n ,Z n ,F n )} C J? 2 (fl) with 

X n G u€ , Z n — F n € i 

If X n — > X strongly and (Z n ,F n ) — (i?, F) weakly in Jzf 2 (r2), then 

Z-FeJtTjt, (4.8) 
Z < F if F n — > F strongly in ^f 1 (fi) and Z n -< F n . (4.9) 

Analogously, for any T > consider {(X n , Z n , F n )} C ,5f 2 ((0, T), Jz? 2 (ft)) with 

X n {t) G JT, Z n (t) - F„(t) G ^ n(t) , Z n (t) -< F n {t) for a.e. t G (0,T). 
7/X„ — > X strongly and (Z n ,F n ) — - (Z, F) weakly in J^ 2 ((0, T), «5f 2 (ft)) 7 then 

Z(t) - F(t) G J# (t) a.e., 
Z(t) -< F(i) a.e. i/F„ — > F strongly in jSf 1 ((0,T), JSf x (n)), Z n (t) -< F n (t) a.e. 

Proof. By assumption, we know that 

Z n {m)<p(X n (m)) dm= [ F n {m)y{X n (m)) dm (4.10) 
for every </? G Cf,(K). Passing to the limit in (|4.10l) we get 

Z(m)(p(X(m)) dm = / F(m)<f(X(m)) dm, 



which yields (|4.8|l since the set o X : <p G C&(R)} is dense in J#x- 
In order to prove f|4.9[) we pass to the limit in the inequality 



tp(Z n (m)) dm < / ^(F n (m))dm 
! in 

for arbitrary convex functions ^; : M — >• R with linear growth, noticing that 

ip(Z(m)) dm ^ liminf / ^!(-Z n (m))dm, 

ip(F(m))dm= lim / ijj(F n (m)) dm. 



n— voo 



(4.11) 



The corresponding inequality for convex functions V> with arbitrary growth at infinity can be obtained from 
(|4.11[) by monotone approximation. 

The time-dependent result follows by applying Ioffe's Theorem. □ 

Theorem 4.4 (Stability of Generalized Lagrangian Solutions). Suppose that F: J€ — > Jz? 2 (£l) is pointwise 
linearly bounded and uniformly continuous. Consider a sequence {X n } of weak Lagrangian solutions with 
initial data 

X n G X and V„ G 

that converges strongly in Jz? 2 (f2) to X £ and V G J%x- Then there exists a subsequence (still denoted by 
{X n }) with the following properties: 

(1) We have X n (t) — > X(t) in _Sf 2 (fi) uniformly on compact time intervals. 

(2) For anyT>0 we have V n — > V in Jz? 2 ((0, T), Jz? 2 (f2)). 

(3) The limit function X is a generalized Lagrangian solution. 
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Proof. Since (X n , V n ) — > (X, V) strongly in J<? 2 (f2) we can find a convex function ip satisfying (12.33)) and 
lirrv^oo ip(r)/r 2 = oo such that 

[X n ] + sC C for all n. 

Here ^ denotes the functional (|2.36[) induced by ip. By Lemma \2.9\ it is not restrictive to assume that ip 
satisfies (|2.34[) . The estimates of Lemma T3. 141 (with W := X) and Gronwall lemma yields 

V[X n (t)] + V[V n (t)] < C T for all t £ [0, T] and all n. (4.12) 

By Lemma [2.101 it then follows that the X n take values in a fixed compact subset of .if 2 (f2) and X n are 
uniformly Lipschitz continuous in J£ 2 {$X). Recall that pointwise linearly bounded operators F are also 
bounded. We can then apply Ascoli-Arzela theorem to obtain a convergent subsequence, which we still 
denote by for simplicity. The convergence is uniform in each compact time interval and the limit 

function X satisfies the same Lipschitz bound. 

Consider now the sequence {Z n } of functions given by Definition 14.11 Since Z n (t) -< F[X n (t)} for a.e. 
t and since F is bounded, (I4.12[) implies that the Z n are uniformly bounded in Jz?°°((0, T), J?? 2 (f2)) for all 
T > 0. Extracting another subsequence if necessary, we may therefore assume that 

Z n ^Z weak* in Jz?°°((0, T), J^ 2 (ft)). 

On the other hand, by uniform continuity of F we have that 

F n := F[X n ] — ► F[X] =: F strongly in if 2 ((0, T), if 2 (r!)). 

Lemma FOl then shows that Z satisfies (14. 2[) . 

The uniform bound on Z„ implies that the maps 

F„(t) := V„ + [ Z n (s) ds for all i > 
J o 

are uniformly Lipschitz continuous in each time interval [0,T] with values in «5? 2 (57). Starting from (|4.1I) 
and applying standard stability results for differential inclusions (cf. Theorem 3.4 in [8], here the strong 
convergence of X n is crucial) , we obtain that X solves 

X(t) +dIjt(X(t)) 3 V+ f Z{s)ds for a.e. i > 0. (4.13) 



o 



In particular, the map X is right-diffcrcntiable in «5f 2 (f2) for each t 0, with right-continuous right-derivative 
V; see Propositions 3.3 and 3.4 in [S]. Therefore gTg|) holds for all f > if is replaced by V(t). We 
may also assume that 

V n — >■ weak* in ^f°°((0, T), jSf 2 (£l)) 
for all T > (extracting another subsequence if necessary). To show that V n — > V strongly in «if 2 ((0, T), ^f 2 (51)), 
we multiply the differential inclusion by 

^((T-t)X„(i)) = {T - t)V n {t) - X n {t) 

and integrate in time over (0, T) x il. Now notice that since X n (t), V n (i) € J%x n (t) an d since dlje(X) C Jtfx 
for all X £ J^, the subdifferential terms vanish after integration over Q. Integrating by parts in the force 
term, we obtain 

J (T-t)\\V n (t)\\% 2{n) dt- (J v n (t,m)X n {t,m)6mjdt 

= -T f V n (m)X n (m)dm - [ (T-t)( [ Z n (t, m)X n (t, m) dm] dt. 
Jn J a \Jn J 

A similar identity holds in the limit. Since the sequence {A„} converges strongly and the sequence {(V n , Z n )} 
converges weakly, we can pass to the limit and get 

lim f (T-t)\\V n (t)\\% 2{n) dt= [ (T-t)\\V(t)\\% 2{n) dt 

for every T > 0. This, together with (|4.12j) yields the desired strong convergence. Therefore there exists an 
/^-negligible set N C [0,oo) such that (up to extraction of a subsequence if necessary) V n (t) — > V(t) in 
^f 2 (57) for every t e [0, oo) \ N. We can then pass to the limit in (|4.3[) written for (X n , V n ) and obtain the 
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corresponding inclusion for (X,V) in (t\,oo) for all t\ G [0, oo) \ N. Since V is right-continuous, formula 
(|Q|) eventually holds for all t x ^ 0. 

(|4.4|) follows by the same argument, first passing to the limit for t\ <E [0, oo) \ N and recalling that by 
(|2U) if K n = Pjt{H n ) and K n — > K, H n — ^ H in j£f 2 (Q) then if = Pj%r(H). □ 

We conclude this section with the main existence result for generalized Lagrangian solutions. As for sticky 
evolutions, its proof relies on the discrete particle approach we will study in the next section, see Remark 
1531 

Theorem 4.5 (Existence of generalized Lagrangian solutions). Let us assume that the force functional 
F: Jff — > Jif 2 (n) is pointwise linearly bounded and uniformly continuous. Then for every couple X E 
and V € Jf?x there exists a generalized Lagrangian solution with initial data {X , V). 



5. Dynamics of Discrete Particles 

We discussed in the Introduction that the conservation law (jl.ip formally admits particular solutions for 
which the density consists of finite linear combinations of Dirac measures; see (|1.6[) above. In this section, 
we will reformulate these solutions in the Lagrangian framework and will prove their global existence. In 
fact, they are Lagrangian solutions in the sense of Definitions 13.41 and 14.11 

For every N € N let us introduce the convex sets 

JV 

W N 



jm e M W : mi > and ^2 mi = lj, 

i=l 

= {a: e 1^: xi x 2 < . . . < iffj. 



For all times t ^ 0, a discrete solution to (jl.ip of the form (|1.6p is therefore determined by a unique number 
N E N and a vector (m, x, v) G x K N x Mr . To find a Lagrangian representation of (|1.6p we consider 
a partition of f2 given by 

i 

=: wo < w\ < . . . < wn '■= 1 where Wi := (5-1) 

2=1 

for i = 1, . . . , iV — 1. Writing Wi := to,) we define functions 

JV W 

X:=^a;ilw 4 and V := ^ttilw,, (5.2) 
i=l i=l 

the (finite dimensional) Hilbert space 

JV 

,Kn ■= {x = J2 x i 1 m : x = g R^} c J? 2 (tt) (5.3) 

i=l 

and its closed convex cone 

N 



\x = 5^3*1^: x = (xi,-- - ,xjy) e K w } cJf c Jz* 2 (ft) (5.4) 



i=l 



Then clearly X G C and V £ , and we easily have 

N N 



Q 



= X#m = ^2 m i ) V = woX, (g«) = y^m^i (5.5) 
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5.1. Discrete Lagrangian solutions. We can reproduce at the discrete level the same approach we fol- 
lowed in Section [31 we can introduce the projected forces 

N 

F m [X] := Pjr m (F[X}) = Va™,*!^, a m<i = f F[X(t)]{m) dm, (5.6) 

which satisfies the analogous of (|3.2p 

j il>(x)f[g](dx) = J ip(X)F m dm, if X G g = X # m as in (JSSJ), (5.7) 

and we can simply solve the differential inclusion 

X(t) + dI jern (X) = V+ f F m [X(s)]ds, X(0) = X 

Jo 



(5.8) 



for given initial data (X,V) G x ^x- Introducing F(t) := V + J Q F m [X(s)]ds = D i=1 !/i(i)lWji we 
end up with the system 

V „ fort)0, (I(0),y(0)) = (X,n (5.9) 

t r (t) = -£< m [A (t)J, 

which is equivalent to (|1.26l) . 

If, e.g., F is Lipschitz, then F m : ,J^ m — > fflm is also Lipschitz and the analogous statements of Theorems 
13.51 and l3~6l hold at this discrete level. In particular, as in (|3.9[) . we have 

V(t) = ^X(t) = P Tx(t) ^ m (Y(t)); (5.10) 

the discrete analog of Lemma T2.4I thus justifies condition (|1.15j) we introduced in the simplified situation of 
a collision of two particles. 

Let us now consider a sequence X n of discrete Lagrangian solutions of (|5.8[) corresponding to initial data 
(X n , V n ) G .%m n x ^m n strongly converging to (X, V) G x Jz? 2 (Sl). We want to show that X n — > X 
locally uniformly in C([0, oo); Jz? 2 (f2)) where X is the Lagrangian solution associated to (X,V). To make 
the analysis simpler, we will assume that the distributions of masses m n give raise by (|5.1[) to suffiently fine 
partitions of the interval (0, 1), i.e. 

for every K G Jlf there exist K n G J^m n such that K n — > K in «5? 2 (£7). (5-11) 

Since J^m n C ^ , (|5.1ip is equivalent to say that the sequence J^m^ Mosco-converge to J?T in the Hilbert 
space Jzf 2 (f2) [TJ Section 3.3.2]. By first approximating C 1 ([0, 1]) functions (which belong to Jlf — J(T) and 
then applying a density argument, it is not difficult to show that (|5.1ip implies a similar property for the 
closed subspaces Ji? mn m «£? 2 (S1), i.e. 

for every H G Jz? 2 (fi) there exist H n G J%n n such that H n — > H in «5f 2 (fi). (5.12) 

Both (|5~TT|) and ([5TT2"]) surely holds if, e.g., 

lim ||m„|| 00 = 0, 

n— >oo 

where for a generic m = (mi, • • • , m^r) G M ff we set ||?n||oo = sup^ 
It is not surprising that we have the following approximation result: 

Theorem 5.1 (Convergence of discrete Lagrangian solutions). Let F : — > Jz? 2 (f2) be Lipschitz and point- 
wise linearly bounded, and let m n G M"" be a sequence satisfying (|5.1ip and let X n G Lip loc ([0, oo); Jtf mn ) of 
discrete Lagrangian solutions corresponding to the initial data (X n ,V n ) G J^m n x ^rn„ strongly converging 
to (X ,V) G Jff x J?? 2 (f2). Then X n — > X locally uniformly in C([0, oo); Jz? 2 (f2)) where X is the unique 
Lagrangian solution starting from (X,V). 

Proof. We cannot directly apply the stability estimates of Theorem 13.61 since the discrete Lagrangian solu- 
tions are associated to convex sets depending on n, so we combine the compactness argument of the 
proof of Theorem 14.41 and a classical stability result for differential inclusion 1, Theorem 3.74] generated by 
a Mosco-converging sequence of convex sets. 
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In fact, we can choose a convex and superquadratic functional ip satisfying (|2.33[) such that 

9[X n ] + < C; 

the estimates of Lemma 13.141 (which can be extended to the discrete case) yield 

V[X n (t)] + #[K(*)] < C T for all t 6 [0, T] and all n. 



Arguing as in the proof of Theorem 14.41 we can find a subsequence (still denoted by X n ) locally uniformly 
converging to a limit X G Lip loc ([0, co); Jz? 2 {Q)) which takes its value in JV. We easily get that F mn [X n ] — > 
F[X] in Lf oc ([0, oo); Jz^ 2 (n)) since for every time t > 

||-Prn„ [X n ] ~ F[X] ||_jf2(f2) ^ C|| A n ~ -^||j^2(fi) + ||-^n. n [X] ~ F[X] 

and F mn [X] = P# mn (F[X]) — ► F[X] in J?\n) by f£T2). 
It follows that 

Y n (t) =V n + [ F mn [X n (s)] ds — > F(t) := V+ [ F[X(s)] ds 
Jo Jo 

locally uniformly in C([0, oo);J>f 2 (f2)). We can then apply (TJ Theorem 3.74] to show that the limit X also 
satisfies the differential inclusion 

X + dIje(X) 3 Y 

and therefore it is a Lagrangian solution associated to (X, V). Since the limit is uniquely determined (by 
Theorem 13. 6p we conclude that the whole sequence X n converges to X. □ 



5.2. A sticky evolution dynamic for discrete particles. In this section we will describe a different 
discrete procedure to construct evolution of a finite number of particles. In the general case, this approach 
will lead to generalized Lagrangian solutions; when F is sticking, we will obtain a sticky evolution which in 
in fact will coincide with the construction we considered in the previous section. 

We already explained the basic idea in the introduction: at the discrete level, a collision between two or 
more particles at some time t 1 corresponds to the impact of the vector x with the boundary dK. N (equivalently, 
of the Lagrangian parametrization X with the boundary of ,%f m in Jtifm)'- in this case, we relabel the particles 
and consider the evolution for t ^ if in a reduced convex cone attached to the new configuration up to the 
next collision. 

In order to get a precise description of the evolution, let us observe that the boundary dK N of the 
cone in R N consists of vectors whose components are not all distinct. For any x e dK N we define 
Ii := {fc: Xk = x{\ for all i = 1, . . . , N. Then there exists a interger A' < A and an increasing map 

a: {I,--- ,A'} —►{!,•■• , A} 

with the property that a(j) = xainl^ij-t for all j — 1, . . . , A'. We set 



and m'jVj := niiVi (5.13) 



for all j and obtain a new state vector (m' ,x' ,v') £ M N x K N x 1^ . In terms of the corresponding 
functions X' £ J(f and V G J4?x' (defined as in ((573])), this means that X' e $f m ' , V S J#x' C J4°m' and 

X' = X, V = Pje xl (V). (5.14) 

Starting from this remark, we can now introduce the precise evolution algorithm for the Lagrangian parametriza- 
tion X. Assume without loss of generality that X does not belong to the boundary of in J^m- We 
construct a map t M> X(t) S J?T as follows: On the time interval [to,ti), where to := and t\ > is to be 
determined later so that X(t) does not touch the boundary of dJffrn in [ta,t{), we obtain functions 

N N 

X(t,-)=J2 x i( t ) 1 W i and V(t,-)=5~Jw i (t)l Wi 
i=l »=i 

by solving the system 

X(t) = V(t), V(t) = P^ m (F[X(t)]). (5.15) 
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Since = J%x(t) i R [A) 1^1)1 we notice that the projection onto Jfffn returns a function that is piecewise 
constant on the same partition on which (X, V) is constant. More precisely, we find 

N 

P*?rn(F[X(t)])=y2 ai (t)l W( where a t {t) := f F[X{t)](m)Am 

for i = 1, . . . , N. Hence (|5 . 1 5[) is equivalent to the system 

x(i) = v(t), v(t) = a(t) for alH e [t ,*i), (5-16) 

which is well-defined. The time t\ is taken as the smallest t > for which X(t) hits the boundary of J^fh 
in M'rn- As explained above, at time t\ we can find an integer N' < N and compute a new state vector 
(rh', st', v') G Jt N ' x X N ' x E^' by ([57TB]) . On the interval i 2 ), with t 2 > ii to be determined, we obtain 

AT' W 

=5^4(<)lw; and V(t,-) = ^,v' j (t)l w , 
3=1 3=1 

by solving (15.151) and (|5.16[) with rh replaced by rh' , the initial condition (a/,u')(£i) := (x',i>'), and the 
new subdivision Wj :— [w'j_ 1 ,w'j) defined by 

3 

=: w' < w[ < . . . < w' N , := 1 where w'j := m' fc for all j = 1, . . . , AT'. 

fe=l 

Again the problem reduces to solving a finite dimensional ordinary differential equation and the time £2 is 
taken to be the smallest t > t\ for which X(t) is in the boundary of Jff N . Then we continue in the same 
fashion. 

We obtain an integer K G N, a sequence of "collision times" 

=: to < ti < . . . < <a--i < tx ■= 00, 

and a pair of functions (X, V) such that 

^x(t) = ^x(t k ^), X(t) = V(t), V(t) = P* xw (F[X(t)]) (5.17) 

for all t G [tk-i,tk) and k = 1, . . . ,K. At collision times the space J#x(t k ) is strictly smaller than Ji?x(t) f° r 
all t < tk, which implies that K ^ N. We have 

X(tk+)=X(t k -), V(t k +) = P^ th) (V(tk-)). (5.18) 

It is easy to check that the monotonicity condition (|3.23|) is satisfied. 

5.3. Sticky and generalized Lagrangian solutions for discrete particles. The next Theorem shows 
that by the algorithm described in the previous section we will obtain a generalized Lagrangian solution 
in the original cone starting from the discrete data (X,V); when F is sticking, this coincides with the 
unique sticky Lagrangian solution. 

Theorem 5.2 (Generalized and sticky Lagrangian solutions for discrete particles). 

Suppose that F : J€ — > Jz? 2 (£l) is uniformly continuous. Consider functions (X, V) of the form (|5.2j) for 
some N e N and (rh, x, v) G M N x K N x R N . 

(1) The curve (X, V) described by the previous section is a generalized Lagrangian solution to (|3.ip with 
initial data (X, V). 

(2) If F is sticking, then (X,V) is a sticky Lagrangian solution. 

Proof. Let us first prove that the map t n> X(t) is a generalized Lagrangian solution with respect to the 
choice 

Z{t) := P^ x{t) (F[X(t)}) for all * ^ 0. 

The fact that V is the right-derivative of X follows immediately from the construction. To prove (|4.3[) it is 
not restrictive to assume t% = 0. We argue by induction on the collision times. In the first interval [to,ti) 
inclusion (14.31) is satisfied by taking the null selection in the subdifferential d!jf(X(t)). 
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Assume now that (|4.3[) is satisfied in for some k. Then 

X(t) = V(t k +) + f Pjr (F[X(s)])ds 

„t 

= (V(t k +) - V(t k -)\ + V(t k -) + / P^ x(s} (F[X(s)}) ds (5.19) 

for any t 6 [t k ,t k+ i), by (|5.1T[) . By induction assumption, we have that 

V(t k -) + £ = V + f " P^ x(s) (F[X(s)]) ds (5.20) 
Jo 

for some £ G dI.x{X{t k )). Combining (|5~T9l> and (|530|) . we obtain 

X(t) + £ + (V(t k -) - V(t k +)) = V + £ P^ x(s) (F[X(s)]) ds 
Because of (|5.17[) . we have that 

V(t k -)= lim h- l (x{t k )-X{t k -h)\. 

h— >-0+ V / 

Using (|5.18[) . we then obtain 

V(t k -) - V(t k +) 

= V(t k -)-P^ x(tk) (V(t k -)) 

= Van h- 1 (x(t k ) - X(t k -h)- P^ x(tk) (X{t k ) - X(t k - h))) 
= ^lim h- 1 (Pj? x(tk) (X(t k - h)) - X(t k - hj) . 

We now use Lemmas 12.41 and 12.61 and conclude that V(t k — ) — V(t k +) G dI,ytf(X(t k )), noticing that NxJf — 
dljc(X) for all X G J£^. Property (|3.23l) implies the monotonicity of the subdifferentials, which are closed 
convex cones. This yields 

e + (v(t k -) - V(t k +)) g dIjr{X{t)) 

for all t e [tfc, tfc + i). Identities (14. 4[) and (|4.6p can be proved as in Proposition 13.81 We conclude that X is a 
generalized Lagrangian solution. 

It remains to show that if F is sticking, then (|3 . 1 1) holds. Because of (|3.23[) . we have that X(t) G J%x(s) 
for all s ^ t. Then Definition 13.31 yields 

j (f[X{s)] - P^ x(s) (F[X(s)]j) ds G Wjr^W). (5.21) 

Adding (|5.21[) to either side of 

V(t) + dI x {X(fy) 3 V + [ Pj? XM (F[X(s)}) ds, 

Jo 

we obtain (|3 . 1 1) . Therefore A is a sticky Lagrangian solution. □ 

We already know that any (even generalized) Lagrangian solution induces a solution of the conservation 
law (jl.lj) . Since for each time t ^ the transport map X(t, •) is piecewise constant, it is easy to check that 
the corresponding solution is in fact a discrete particle solution: the density /momentum is of the form (|1.6[) . 

Remark 5.3. Notice that piecewise constant functions as in (|5.2p are dense in Jzf 2 (f2), so we can approximate 
any given initial data and then combine the existence result 15.21 with the stability Theorem 14.41 to get the 
proof of Theorem 14.51 

Remark 5.4. The proof of Theorem 13.111 follows by a similar approximation argument. By Theorem 13. lOl it 
is sufficient to show that any Lagrangian solution A with A G J(f and V 6 J%x satisfies 
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That is, if X is constant on some interval (a, /3) C f2, then X(t) remains constant on (a,/3) for all times 
t ^ 0. We approximate (X, F) by a sequence (X n , V"„) of the form (|5.2[) such that X n is constant on (a, (3). 
Since this property is preserved by the discrete Lagrangian solution constructed in Theorem l5.21 the stability 
estimates of Theorem 13.61 show that the limit function X(t) is still constant on (a, 0). 



6. Global Existence in Eulerian coordinates 
Theorems 13.61 13.111 13.131 and 14.51 of the previous sections 

immediately translate into global existence results for the Euler system of conservation laws (|l.lj) . Before 
stating some of the related results, let us explore in more detail the relation between the force functionals 
f[g] in (jl.ip and their reformulation in the Lagrangian framework. 

6.1. The Eulerian description of the force field. Let us first introduce the space 

,% (R) :={( ft!) ) :(? e^(l),^ Jz? 2 (R, £>) j. 

For all {g l ,v i ) <E -%(R) with i = 1 ... 2, we then define [10 §2] 

D2\Xq\,Vx), {g 2 ,v 2 )j := max|w / 2 (£'i I Q2), U 2 [{qi,v 1 ), (02,^2))}, 

where W 2 is the Wasserstein distance and U 2 denotes the semi-distance 

Ui((gi,Vi), {g 2 ,v 2 U := / \vi(x) - v 2 (y)\ 2 g(dx, dy) 

Wl (X ei (m))-« 2 (X e2 (m))| 2 dm. (6.1) 

n 

Here g £ r opt (p 1 , g 2 ) is the unique optimal transport map between the measures g\ and g 2 . It can be 
expressed in terms of the transport maps defined in (|2.4p ; see (|2.6p . The sequence {{g n ,v n )} converges to 
(g, v) in the metric space (^(R), D 2 ) if and only if W 2 (g n , g) — > 0, if g n v n — ^ gv weak* in ^(R), and if 

\v n \ 2 Q n — > J M 2 g- 

We refer the reader to [I7J Prop. 2.1] and to [3] for further details (see in particular Definition 5.4.3). 

We consider a continuous map (with respect to the Wasserstein topology in 3* > 2 (R) and the weak* topology 
on ^f(R) induced by C 6 (R)) 

/: <3» 2 (R) — ^^(R), f[g] = f e Q, f e eJ? 2 (R,g), (6.2) 

with the property that f[g] is absolutely continuous with respect to g £ ,^ 2 (R): fg is the Radon-Nikodym- 
derivative of f[g] with respect to g and assume that f B £ „Sf 2 (R, g). 

Definition 6.1 (Boundedness). We say that a map /: 3 g 2 {K.) — > .-#(R) as in (j6.2[) is bounded if there 
exists a constant C ^ such that 

ll/Jl«(R lC ) < c(i + / M 2 dg) for all e £ 5» 2 (R). 

We say that / is pointwise linearly bounded if there exists aC p )0 such that 

\f e (x)\ < Cp^l + |x| + y |x| dg) for a.e. x £ R and all g £ ^» 2 (R). 

Definition 6.2 (Uniform continuity I). We say that a map /: ^ 2 (R) — > ^#(R) as in (|6.2I) is uniformly 
continuous if there exists a modulus of continuity lu such that 

U 2 [(giJ ei ),(g 2 J g2 )) ^Lu(w 2 (g u g 2 )^ for all gi,g 2 £ ^2 (M) . (6.3) 

In the case w(r) = Lr for some constant L and all r 0, we say that / is Lipschitz continuous. 
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As discussed in Section 12. 1[ there is a one-to-one correspondence between measures g G 3^2 (K) and 
optimal transport maps X E ,5f 2 (£2), given by 

X EJf and X # m = g. (6.4) 

We now want to construct a functional F; Jif — > Jff 2 (fl) such that 

/ (p(x) f[g](dx) = / <p(X(m))F[X](m)dm for all <p e C 6 (R), (6.5) 

whenever (X, g) are related by (|6 .4[) . One possible choice is to set 

F[X] := f g o X for all (X, g) satisfying (jUl), (6.6) 

which easily gives F[X] G «5f 2 (f2). Then the boundedness and continuity assumptions on the functional / 
in Definitions 16.11 and 16.21 translate immediately into the corresponding properties for F in Definitions 13.11 
and !3.2l It can be useful, however, to also consider different choices for F. 

Definition 6.3 ( Uniform continuity II). We say that a map /: ,^ 2 (R) — > ^#(M) as in (|6.2p is densely 
uniformly continuous if (I6.3[) holds for measures that are absolutely continuous with respect to C 1 with 
bounded densities. We define dense Lipschitz continuity similarly. 

Lemma 6.4. ///: ,^ 2 (R) — ► ^(R) is densely uniformly continuous, then there exists a unique uniformly 
continuous map F: X — > jSf 2 (Q) such that holds for all (X, g) satisfying ([BH) . 

Note that ((63]) imlies that f e o X = ?^ X {F[X\) for all (X, g) with ([Of . 

Proof. We denote by J^ C g the dense subset of whose elements are C 1 (f2)-maps with strictly (hence 
uniformly) positive derivatives. For every X £ J£^ e g the push-forward g := X#m is absolutely continuous 
with respect to the Lebesgue measure C 1 and has a bounded density. We can then define 

F[X] := f e o X for all X e -X cg . (6.7) 

Applying definition (|6.1I) and (|6.3I) we obtain 

||F[Xx] - F[X 2 ][[^2 (fi) s: -X 2 ||^ 2(n )) for all X l5 X 2 e ^ cg . 

Then F can be extended to all of Jtif by density. One can check that this functional satisfies (|6.7p , therefore 
it is uniquely determined by /. □ 

Definition 6.5 (Sticking). Let /: ^ 2 (R) — > ^#(R) be densely uniformly continuous and let F be the 
functional from Lemma 16.41 We say that / is sticking if F is sticking. 

6.2. Existence results and examples. We state here a simple example of possible applications of the 
previous Lagrangian results; for the sake of simplicity, we omit to detail all the information which could 
be derived by the finer structure properties and by the a priori estimates we obtained for the Lagrangian 
formulation. It is worth noticing that all the solutions can be obtained as a suitable limit of discrete particle 
evolutions. 

The first statement follows by Theorem 14.51 the second one by Theorem 13.61 Theorems 13.101 and 13.111 
yields the last assertion. 

Theorem 6.6 (Global Existence). Let us fix g e ^* 2 (R) and v £ «5f 2 (R, g).. 

(1) Suppose that the force functional f : ^ 2 (R) — > J?(R) is pointwise linearly bounded and densely 
uniformly continuous. Then there exists a solution (g,v) of the conservation law with initial 
data (g, v). 

(2) If f : ^ 2 (R) — > ^(R) is densely Lipschitz continuous, then there exists a stable selection of a 
solution (g,v) of (|1.1[) with respect to the initial data (g,v) in (5 2 (R),D 2 ). 

(3) If f : ^ 2 (R) — > ^(R) is densely Lipschitz continuous, and sticking, then there exists a stable sticky 
solution (g,v) of with initial data (g,v). The map St ■ {g,v) i— > g(t, ■), v(t, •)) is a semigroup 
in (5%(R), D 2 ). 

We finish the paper by giving a number of examples of force functionals. 



STICKY PARTICLE DYNAMICS WITH INTERACTIONS 



.37 



(6.8) 



Example 6.7. Let v: K — > K be a continuous function satisfying 

\v(x)\ s$ C„(l + \x\) forallxgR, 

with C v some constant. Then the operator defined by 

f[g]:= gv for all g g ^ 2 (R) 

is pointwise linearly bounded; it is Lipschitz continuous if v is a Lipschitz function. Note that f[g] is the 
Wasserstein differential of the potential energy (see [J] ) 



V[g] := / V{x) g(dx) where v = V . 
Jn 

Example 6.8. Let w: R — > R be a continuous function satisfying (|6.8|) . Then 

f[g] := g{w * g) = q(J^{- - y) g(dy)^j for all g g ^ 2 (R) 

is pointwise linearly bounded, since 

|(«;*e)(a;)| < C w / (1 + |x - y\) g(dy) ^ C w ( 1 + |x| + / \y\ g(dy)) 

It is Lipschitz continuous if w is a Lipschitz function. In fact, writing f g :—w-kg for all g G ,^ 2 (R), we have 
that 

|/eiO*0 -/ e2 (y)l = / w(x-x')g 1 (dx')- J w(y-y')Q2(dy') 

(w(x - X 7 ) - w(y - y')) £>(da;', dy') 

^L^\x-y\+J^ \x' -y'\Q(dx',dy')\ 
where L ^ is the Lipschitz constant of w and £ g F opt (f?i, £2)- This implies 

%((f?i,/ ei ),0?2,/ e2 )) <4LW 2 (gi,g2) forall0i,0 2 g^2(R). 
Note that /[g] is the Wasserstein differential of the interaction energy (see [4]) 

W(x — y) g(dx) g(dy) where w = W' . 



Example 6.9. Let us consider the previous example with the Borel function 

f 1 if x > 
w{x) := I if x = , 
[-1 ifx<0 

which corresponds to := \x\. To show that f[g] is continuous, note that 

f e (x) = m g (x) + M e (x) - 1 for all 

where m g (x) := 00, x)) and M s {x) := g{{— 00, x]) as in (|2.3[) above. Up to rescaling and adding 
constants, the function f g is the precise representative of the cumulative distribution function of the measure 
g. For convenience, we define 

f e {x) := f e (x) + 1 for all x e M, f[g] := / e £>. 

We now introduce the sets 

J := [x g R: g({x}) > o} and J" e := |J (m e (x), M g (x)) . 

Note that J e is at most countable. If X g is defined by (|2.4j) . then 

X e (m) = x for all m G [m g (x), M g (x)] and x £ J g , 
f g (X g (m)) = 2m for all m Ett\J g . 
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For any ip G C(,(R) we have 

<p(x) f[ Q ](dx) = I <p(X e (m))f B (X B (m))dm 



M e (x) 

ip(X e (m))f e (X e (m))dm+J2 / <p(X B (m))f e (X e (m)) dm 



_(> 



2 / (p(X g (m))mdm+ V (m%(x) - m 2 g (x)) ip(x) 
2 / <p(Xg(m)) mdm + 2 y I tp(x)mdm 

Jn\J„ r Jrn„(x) 



-M e (x) 

'CI\Jq x£j„ J ' m ei x ) 

= 2 / ip(X g (m))mdm 



(p(x)f e (x) g(dx) = / tp(X g (m)) (2m - 1) dm 



It follows that 



Then the map / is pointwise linearly bounded because \f g (x)\ ^ 1 for all x G R. It is continuous since 

g n — > g in ^2 (R) implies that X Bn — > X e in J2? 2 (il). It is densely Lipschitz continuous since the 
associated functional F is given by 

F[X](m) := 2m - 1 for all m G Q, (6.9) 

which does not even depend onle anymore. 

Example 6.10. For er G Jz?°°(IR) let q g be the solution of (recall (fO]) ) 

-9L9 e = A( e - ( r). (6.10) 

Then g e is locally Lipschitz continuous and its (opposite) derivative a g := —d x q g is locally of bounded 
variation. Choosing its precise representative we then define 

f[e] ■= Qa g for all g G ^ 2 (R)- (6.11) 

Setting (5o-(x) := cr(y) dy it is not difficult to check that 



(x) = -A(-(m e (x) + M g (x) - 1) - Q CT (a;)) 



for all x G 



so that the associated operator F is given by 



F[X](m) = -\{-{2m - 1) - Q CT (X(m))) for all to G O. 

This corresponds to the Euler-Poisson system discussed in the Introduction. For simplicity, let us consider 
consider the case when a vanishes. 

Sticky solutions for the attractive Euler-Poisson system. In the attractive case (when A > 0) the functional 
F[X] is sticking: Let fix be defined by (|1.32[) and let (a, /?) C fix be a maximal interval. Then Pj^ x (F[X]) 
is constant in (a, /3) and equal to its average over the interval. We define 

E(m):=f (F[X](m) - P# x (.F[X])(m)) dm for all m G (a, 0). (6.12) 



Then S(a) = 5(/3) = 0. Since A > 0, 3 is concave and, we obtain that S(m) ^ in (a,/3). By Lemma [ 
we conclude that the functional F is sticking. Sticky Lagrangian solutions to the Euler-Poisson system (jl.ip 
(thus obtained as limit of sticky particly dynamics) are therefore unique and in fact form a semigroup in the 
metric space (.%(R),D 2 ) by Theorems EH [331 and E H 

We can then apply the representation formula (|3 . 26[) and (|2.9|) to obtain the following result: 
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Theorem 6.11 (Representation formula for attractive Euler-Poisson system). The unique sticky Lagrangian 
solution of the Euler-Poisson system (X ^ 0) corresponding to initial data (g, v) with g = X#m, X 6 , 
and V = v o X , can be obtained by the formula 

g(t,-)=X(t,-) # m, X(t,m) = ^-X**(t,m) (6.13) 

am 

where X**(t,m) is the convex envelope (w.r.t. m, see (|2.10l) ) of 

f m / - - t 2 \ 

X(t,m):=J {X{l)+tV(l)-\—(2l-l)}dl (6.14) 

Notice that when A = we find the sticky particle solution of [T7] . 

Lagrangian solutions for the repulsive Euler-Poisson system. In the repulsive case A < the function S 
defined in (|6 . 1 2[) is convex and vanishes at the endpoints of (a, /3), thus S(m) ^ for all m G (a, (3), and the 
map F does not satisfies the sticking condition. In this case (|6 . 1 3[) - ()6. 14[) may be different from the solution 
given by Theorem 13.61 

Here is a simple example for A = — 2: consider the initial condition 

X(m) :=m-l/2, V(m) := - sign(m - 1/2), (6.15) 

for which (|6.14[) yields 

X(t, m) = -(1 + i 2 )(m- 1/2) 2 -t\m- 1/2| - c(t), c(t) := -(1 + t 2 - At). (6.16) 

2 8 

It is easy to check that 

(X(t,m) if \m-l/2\>5(t), t 

x ^ m )=(__ n £ FT _ c(t) if | m -i/ 2 | <«*(*), where 5( ' ):= TT^' (6 ' 17) 

so that X(t, ■) is the piecewise linear continuous map 

, (x(t,m) if [to - 1/21 > 5(t), , -. , 9 , 
X(t,m) = < v ; , , where Xft.m) := — #(t,m). 

v ' ; |0 if \m-l/2\<S(t), y ' ' dm v ' ; 

If we eventually introduce 

d d 2 

Y(t, m) := —X(t, m) = Q^^{t, m) = 2t(m - 1/2) - sign(m - 1/2), 
recalling (|3.6[) X is a Lagrangian solution if and only if 

Y(t, •) — X(t, ■) G a/jfprft .)) a.e. in (0,oo). 
By Lemma 12.31 we obtain the equivalent condition 

d 

— (X(t,m)-X**(t,m))>0 in - S(t) < m < S(t), (6.18) 
which is not compatible with (|6.16l) and (I6.17[) : to see this, fix e.g. < 5 < 1/2, := 1/2 + 5, and 

t± ■■= 1±V 28~ 4S2 > 80 that <*(*±) = 6 < for ever y i e (*-)*+)• We thus have 

A'(i ± ,TO*)-A'**(i ± ,TO*) = 0, #(t,m*) - X**(t,m,) > for i_ < t < t+, 
which contradicts (|6.18|l . 

7. Convergence of the Time Discrete Scheme of Section 11.41 



In this section, we establish the convergence of the time discrete scheme of Section 11.41 Since the proof 
does not substantially differ from the one provided in [7] for order-preserving vibrating strings, we only 
sketch the main steps. The key point is the non-expansive property of the time-discrete scheme. Indeed, we 
first observe that the rearrangement operator, even in the periodic case, is non-expansive in J£ 2 (Vl). More 
precisely, we have that 



1 rl 



\Y*(m) - Z*{m)Y dm s; / \Y(m) - Z{m)\ 2 dm 
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for all pairs (Y, Z) of maps such that Y — id and Z — id are 1-periodic and square integrable. Next, we see 
that the harmonic oscillations (|1.39p are isometric in phase space for (X(t,m) — m,V(t,m)), for each fixed 
m. Let {X Ttni V T _ n ) 1 (Y TM ,W T ^ n ) be generated by the time-discrete scheme. Then 

|Xr,n+l - ^r,n+l||j5?2(n) + ||^r,n+l ~ Wt,»+1 || SC 2 (n) 

^ ||^r,n+l _ ^r,n+l|||f2(o) + ||K,ri+l - Wr,n+l||if2(J2) 
= \\X T ,n — ^r,n|||f2(Q) + ||K",n ~ ^r.n\\%2^Qy 

Since (X = id, V = 0) is a trivial solution of the scheme, we immediately get 

H-^r.n+l - id ||^2( n ) + ||K,n+l|||f2( n) ^ \\X - id||| f .2 (n) + ||V"||j^2( n ). 

Because the scheme is translation invariant in m and (discretely) in n, we easily deduce the strong compact- 
ness in Cj(jSf^) of the discrete solutions, linearly interpolated in time, for each 1-periodic initial condition 
(X— id, V), first in Jf? 1 and then (by a density argument, using the non-expansive property of the scheme) in 
Jz? 2 . Let us now examine the consistency of the scheme. To do that, let us compare a solution of the discrete 
scheme to any smooth test function m — > (Y(m), W(m)) where Y is nondecreasing and (Y(m) — to, W(m)) 
is 1-periodic. Since the rearrangement operator is non-expansive and Y = Y* is nondecreasing, we first get 

\\X T , n+ i-Y\\% 2(m + IIK n+ i-W» 2 



1 1 {X T , n (m) — to) cos(t) + V T , n (m) sin(r) — (Y(m) — m) | 2 
+ \ (X T _ n (m) — to) sin(r) — V T . n (m) cos(t) + T / F(to)| 2 | dm. 



One can then check that 



\X T , n+1 - Y\\^ 2 ^ + ||K,n+l - W / Illf2( n ) 

*S \\X T>n -Y\\% 2{n) + \\V T , n -W\\% Hn) 



+ 2t f {(X T . n {m) - Y(m))V T , n (m) - (X T . n (m) - m)(V T , n (m) - W(m))} dm + 



KT 2 , 



with constant k depending only on the test functions (Y,W) and the initial data (X,V). Clearly, this 
estimate is consistent with the differential inequality 

±{\\X(t, •) - Y\\% 2{n) + \\V(t, •) - W\\% Hn) } (7.1) 



^2 ' 



{(X(t,m)-Y(m))V(t,m) - (X(t,m) - m)(V(t,m) ~ W{mj)\dm, 



valid for all pair of 1-periodic functions of form m n> (Y(m) — m, W(m)) with Y nondecreasing, which is 
nothing but the "metric formulation" of (|1.35l) . Indeed, for a.e. t ^ fixed, by choosing Y — X(t, •) and 
W = V(t, -)±Z for arbitrary 1-periodic Z € i? 2 (ft), we find that 

V(t,m) + X(t,m) - m = 0; 

cf. (TE3SJ). On the other hand, by choosing W = V(t, •) and Y = rcsp. Y = 2X(t, •), we obtain 

l 

X(t, m)(X(t, m) - V(t, to)) dm = 

/ Y(m)(X(t, to) — V(t, m)) dm for all Y nondecreasing with Y(m) — to 1-periodic. 
Jo 

This implies precisely that —X(t, ■) + V(t, •) € dIj%r(X(t, •)), which gives (| 1 .35|) . This concludes the proof of 
convergence for the time-discrete scheme. 
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